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ON THE INFLUENCE OF THE ANISOTROPY OF THE BASE ON 
THE RESISTANCE OF A ROLLING CYLINDER 


Bv Japwıca HALAUNBRENNER 
Department of Physics, Cracow Polytechnical College 
(Received June 8, 1957) 


In this paper the directional dependence of the resistance of a cylinder rolling on 
anisotropic surfaces is determined experimentally for the case of artificial anisotropy 
produced by machine working and of natural anisotropy of wood and NaCl monocrystals. 


Experiment shows that a cylinder rolling on a horizontal plane diminishes its 
speed and comes to rest after some time. This can be described phenomenologically 
by introducing the notion of rolling “friction“, or still better, of rolling resistance, 
as a force directed oppositely to the velocity of the cylinder axis. It was already shown 
by Coulomb that this resistance T is directly proportional to the force N pressing the 
cylinder to the base, and inversely proportional to the cylinder radius R. 


N 


The proportionality coefficient k, called the coefficient of rolling friction, has the 
dimension of length; its value depends on the material of both bodies, the state of 
their surfaces, the rolling speed, and, in the case of anisotropic bodies, on the direction 
in which the cylinder is rolling. 

The subject of this paper is the experimental determination of the dependence 
of the rolling resistance on the direction of motion of the cylinder in the case of an 
anisotropic base. We distinguish here the anisotropy of the surface itself resulting 
from the application of cutting processes to an isotropic material, and anisotropy 
appearing throughout the entire mass of the base, such as occurs е. g. ш wood and 


in monocrystals. | 
In order to choose properly the experim 
which physical phenomena contribute to the loss o 


der. We list below the processes involved. E : 
1. Compression of Ше surface asperities and jumps of cylinder over them. 


2. Consolidation (strengthening) of the material of the cylinder and base. 
j (83) 
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ental conditions it is of use to consider _ 
f kinetic energy of the rolling cylin- 
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3. Abrasion and wear of the material 
4. Cohesion of the cylinder and the base: Tomlinson (1929) in his molecular 
theory of friction sees as the essence of rolling friction precisely the cohesive forces 
between the cylinder and the base. 
5. Sliding friction; the rolling cylinder, pressed to the base, becomes deformed 
and deforms the base. The parts of the cylinder in the area of contact are compressed, 
and parts of the base are under tension. As 
the cylinder rolls the area of contact moves, 
the parts of the base not lying on the central 
line and released from compressive stress 
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slide over the surface of the cylinder; the 
parts of the base surface coming under 
the cylinder similarly slide over the sur- 
4 face of the cylinder. Reynolds (1875) and 
2, Föppl (1947) ascribe a decisive role in 
the dissipation of energy to sliding friction 


i 


Fig. 1. Velocity of the points of the base in and calculated the relation between rolling 
the area of contact with the cylinder resistance, and the coefficient of sliding 
friction. 

6. Elastic deformation of the base. Although at first glance it might seem that 
in the case of a perfectly elastic base the work done in compression by the front part 
of the cylinder is given up during the straightening out under the rear side, after 
considering, however, the velocities of the points of the cylinder in the area of contact 
(Fig. 1), we see that the velocity distribution, and therefore the pressure of the rolling 
cylinder on the base, is not symmetrical with respect to OS as in the case of a cylinder 
at rest. The rear part of the cylinder “moves away“ from the base as the latter straight- 
ens out, the front part experiences inertial resistance from the base in addition to 
elastic resistance. This results from the moment of forces hindering the rotation. 
Another factor in the dissipation of energy is the imperfect elasticity of the cylinder 
and the base owing to the internal friction of solid bodies; even in the region of 
proportional deformations a body free from external forces does not give up all the 
work done on it. 


7. The production of elastic waves in the base. These waves carry away consider- 
able energy. Thus e. g. the vibration of rails is observed at a distance of several km 
in front of and behind a moving train. 

8. In the case of a large stress in the metal a plastic deformation takes place. 


9. Finally, the rolling cylinder loses energy in overcoming the resistance of the 


air and the sound waves produced in it. У 


10. In the case of dielectrics, the rolling cylinder loses energy in overcoming the 
different electric polarization of the cylinder and base. 

The percentage of each of these factors in the total resistance depends on the 
type of material, the pressure, the state of the surfaces, and the rolling speed. It seems 
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that in the case of metal, at large pressures, the dominating role is played by plastic 
deformation; in the case of rubber, by internal friction. The phenomenon of cohesion 
can only occur with very clean surfaces. Let us now consider which of the phenomena 
1—10 can affect the directional dependence of the friction coefficient. 


At small pressures, the surface and its asperities will have a dominating in- 
fluence. Abrasion of the material, especially in the case of metals, does not play a big 
role although the abrasive properties of crystals show a marked dependence on the 
direction of grinding (Kuznetsov 1952). At larger pressures the directional differences 
in the elastic coefficients should be felt. At pressures exceeding the yield point it 
should be expected that the directional dependence of the friction disappears, and 
similarly at greater velocities (of the order of 1 m/sec). 


In order to determine the coefficient of rolling friction in the case (1) of artificial 
anisotropy (produced by cutting processes) the inclined plane method was used; 
in the case (II) of natural anisotropy, the rolling vibration method. 


L 


Let us carefully place an isotropic cylinder on an inclined plane with the axis 
perpendicular to the direction of the greatest slope of the incline. We neglect the 
resistance of air for the time being. If there was no friction force T, the cylinder would 
slide down the incline with an acceleration g sin a. If we assume that the reaction of 
the base В (with components М, T) is applied at the point of contact O between the 
cylinder and the base (Fig. 2) and if we assume rolling without sliding, then the ac- 
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Fig. 2. Force diagram of a cylinder rolling on an inclined plane 


celeration of the cylinder axis will have the value: x = 3 g sin M and the friction 
force will be T = $ т g sin 4. The expressions for the acceleration of the axis and 
the friction force do not contain parameters characterizing the material and the surface 
structure of the cylinder and inclined plane. The work done by the force of ae 
equals O, since the velocity of the point of contact is zero. The above description o 


# 
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the rolling cylinder therefore does not correspond with reality. We obtain a better 
description if we consider that both the cylinder and the base undergo deformation 
and are in contact with each other over a certain area and we shift the point of appli- 
cation of the reaction В (T, №) by a distance Е in the direction of motion of the cylinder 


А 


—/ 
GE P 
Fig. 3. Force acting on a cylinder rolling down an inclined plane. The point of application of the 
reaction R is shifted by a distance k in the direction of motion of the cylinder 


(Fig. 3). We denote by r the radius of the cylinder, by J, its moment of inertia with 
respect to the horizontal axis passing through the centre of gravity, by w the angular 
velocity of the cylinder. The equations of motion of the cylinder have the form 


mx = та sin a—T 
my=me cos a—N=0 
J.o=Tr—Nk>0 


Adding the condition of rolling without sliding: œ = xjr and setting k = const. and 
T = const. we have 


ER sin & 2 
SR — 74508,0 2 


The cylinder rolls with a smaller acceleration than before. 
The size of the coefficient k characterizes the material and the state of the surface 


of both bodies. If we determine 2, a, and г experimentally, we can calculate Ё: 


2 g cosa 

In the measurements were used ground steel cylinders 40 mm long and 10—20 mm 
in diameter. The base was in the form of soft steel plates 15 cm x 14 cm in size. The 
plates were worked by a planing machine and placed on the incline, the direction of 
the ridges from the machining operation making an angle of y, with the direction 
of the greatest slope, where 0? < y < 360°. The angle of the incline with respect 
to the horizontal was 2°—5° in order to exclude sliding; the angle was measured with 
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a goniometer to an accuracy of 2’. Under these conditions the cylinder rolling time 
was about 1 sec. A horizontal disc rotating with a constant speed was used to record 
the time. А pen rigidly fixed to the armature of an electromagnet drew a circle on the 
disk. The cylinder, energizing the coil of the electromagnet at the beginning of its 
motion and de-energizing it at the end of its motion, moved the pen several mm. The 
rolling time was calculated to an accuracy of 1% from the length of the recorded arc. - 
A greater accuracy in the measurement of the time was not necessary because the spread 
in the results, as is usually the case in friction experiments, was several times greater. 
The acceleration of the cylinder axis was calculated from the formula 5 = 2. 

This method is not marked by great sensitivity, but has, however, the advantage 
that the rolling cylinder comes in contact with the base over a considerable area, 
which decreases the influence on the result of places not having the average geometry 
and gives a more certain value for the coefficient k. 

The table below contains the results of the observations and calculations for 
a ground steel cylinder 40.2 mm long, 11.2 mm in diameter, and mass of 30.6 rolling 
over an iron base having a smoothness of h = Ви and p = 0.214 mm; the angle 
of the incline with the horizontal is œ = 2?30'. The period of rotation of the disk is 
T= 0,756 sec. 

With the grooves of the base parallel to the motion of the cylinder (у = 0°) | 
4 successive measurements gave for the angle of rotation of the disk during the rol- 


ling of the cylinder the values ф = 470°, 497°, 488°, 477°; hence PES 483° and t = 


— 0.756 sec: = = 1.01 sec. Since S = 14.1 em,. then 
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The smallest value k = 0.0011 cm corresponds to the case in which the cylinder axis 
is perpendicular to the direction of the machining ridges, the largest value 
Е = 0.0044 cm, to the case in which the axis is parallel to the ridges. Rotation 
of the base by 180° has a well resolvable effect. This is explained by the directional 
asymmetry and the shape of the scratches in the furrows left by the planer 
knife. 


Ше 
The case of natural anisotropy was investigated for pine wood and sodium chloride 


crystals. For this purpose there was employed the very sensitive method of a vibrating 
cylinder with a lowered centre of gravity and resting on a horizontal base (Fig. 5). 


Fig. 5. Diagram of the arrangement for optical recording of the vibrations of the cylinder 


The centre of gravity of the cylinder was lowered by placing lead sheets on it 
eccentrically. A very small totally reflecting prism is located on the cylinder 
axis. It directs a beam of light onto a rotating cylinder covered with photo- 
sensitive paper. Lens b gives an image of a distant “point“ source on the surface 
of the cylinder. 

The cylinder displaced from the equilibrium position executes damped vibrations, 
successive amplitudes of which can be measured on the photogram. As a measure of 
the resistance we take e. g. the number of periods during which the amplitude is 
reduced to a value n times (e. g. 2) smaller than some initial value. This method is 
exceedingly sensitive, too sensitive to be used for investigating the rolling resistance 


Resistance of a Rolling Cylinder 89 


in the case of the worked surface. It is seen from the photogram (Fig. 6) how the cylinder 
rolling perpendicularly to the furrows changes its amplitude in jumps and rebounds 
from their sides with a rocking motion of considerably shorter period. 

Figs. Ta and 7b represent photograms of the damping of the vibrations of a polished 
steel cylinder of 11.0 mm diameter loaded with lead disks to a mass of 1.7 kg on a pine 


Fig. 6. Fig. 7a. 
d base for the case a) with the axis parallel and b) with the gi perpendicular 
ES S з . . . = 
er b s of the wood. The difference in damping attains 50 ТА s ~ ledega 
= d P d example of directional dependence of rolling resistance т + è mob 
E i - tal. A па 
t as of the base is the rolling on the surface of a monocrysta 
natural an 
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NaCl monocrystal in the form of a cube with 5 cm edges was used. Since the natural 
faces of the cube have a layer-like structure one of the faces was polished. А vibrating 
steel cylinder of 10 mm diameter was placed first parallel to one edge of the polished 
face and then along its diagonal. Since the expected anisotropy could be caused 
primarily by the difference in the elastic coefficients of the crystal in each direction 
the cylinder was loaded on its axis with lead disks to a total vibrating mass of 4 kg 


Д— = 


Fig. 8a. Fig. 8b. 


in order to increase the elastic deformation. Fig. 8a shows the photogram of the 
vibrations of the cylinder with its axis parallel to the crystal edge; Fig. 8b, with the 
axis making an angle of 45° with the edge. The effect is marked: in the case of the 
cylinder axis parallel to the crystal edge, the amplitude (double) after 30 vibrations 
drops from 33 mm to half the value; in the case where the cylinder laid along the 
diagonal — after 50 vibrations. The effect can be reproduced well at different places 
on the crystal and testifies to the great sensitivity of the method. It permits the discovery 
of elastic anisotropy in an optically isotropic crystal. 
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| КРАТКОЕ СОДЕРЖАНИЕ 


3 Я. Галаунбреннер, Влияние анизотропии основы на сопротивление катящегося 
цилиндра. 

| В этой работе установлено экспериментальным путем зависимость сопротив- 
ления от направления при качении цилиндра по анизотроповым поверхностям, 
в случае искусственной анизотропии, вызванной машинной обработкой и есте- 
ственной анизотропии дерева и монокристаллов NaCl. 
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MICROSCOPIC OBSERVATIONS OF THE SURFACE OF InSb 
MONOCRYSTALS 


Ву J. AULEYTNER AND В. KOŁAKOWSKI 


Institute of Physics of the Polish Academy of Science, Warsaw 
(Received Juni 17, 1957) 


Presented in this paper are the results of microscopic observations of some etched 
and cleaved surfaces of.InSb monocrystals. Surfaces with crystallographic indices (100), 
(111), and others were investigated. It was shown that the shape of the etch pits is different 
for the various crystallographic surfaces. Two types of etched figures were observed; 
they differ in size and in density. Suitable etching agents were chosen to reveal the figures. 


The results are given of the microscopic observations of etched and cleaved 
planes of InSb monocrystals. Planes with crystallographic indices (110) and (111) 
and others were investigated. It was shown that the shape of the etch pits is different 
on the various crystallographic planes. Two types of etched figures differing in- size 
and density were observed. Suitable etching agents were chosen to reveal the figures. 

The purpose of this work was to investigate the surface structure of indium anti- 
monide monocrystals. In the investigations cleavage planes and ground and polished 
surfaces were etched by various methods. The observations were made by means of 
a (Reichert) metallographical microscope. A magnification of from 120 x to 1000 x 
was employed. The orientations of the planes were determined by the Laue method. 
A fine focus X-ray tube with a focal diameter of the order of 30 y (Auleytner et al. 
1957) was used to determine structural defects in the monocrystals. 

` The surfaces to be observed were prepared in the following way: the suitably cut 
surface was ground by alundum, polished with felt, and then immersed into a saturated 
К] solution, after which it was etched for 10 seconds in etching mixture Nos. l or 4 
and for about an hour in etching mixture Nos. 2 or 3. The etching agents used are 


listed at the end of this paper. 
The etched surfaces were rin 


wards dried. The cleavage planes remained untreated. he ii 
In the first part of the work the aim was to carry out an optical investigation of the 


cleavage plans. These investigations are important because they permit us, in part, 
to become acquainted with the kinetics of the cleavage of monocrystals (Gilman 1956). 
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sed in methyl alcohol and then in acetone and after- 
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Indium antimode crystals show very good cleavage along the (110) plane. 

A series of photographs of the planes (111) and (110) were made. 

Figure 1 shows a general view of the observed region of the cleavage (plane (111)). 
The bright band at the centre is a number of parallel lines. By using a greater 
magnification (Fig. 2) the shape is seen to approximate a twisted cord. A number of 
parallel lines is visible on the triangular region to the right. Figs. 3 and 4, representing 
a considerably greater magnification, disclose their finer structure. Figs. 5 and 6 
represent steps analogous to those described by Gilman (1956) characterizing the 
kinetics of crack propagation. Fig. 7 discloses the laminar structure of the steps on 
the surface in which the crack developed. 

The central region in Fig. 8 shows the presence of large stresses in the crystal 
lattice which were produced by irregularities (pits) in the cleavage plane. 

On the (110) cleavage planes are visible lines whose density according to Gilman 
(1956) is inversely proportional to the crack propagation velocity (Figs. 9 and 10). 

The etched surface making an angle of 7” with the plane (111) shows etch pits 
(Figs. 11 and 12). Their mean density is 105/cm?. Their shape is the same as on the 
analogous surfaces of germanium (Auleytner et al. 1957 and other papers). In the 
above, as well as in the other reports, they were regarded as the place of localization 
of single dislocations. In places with larger lattice stresses, and therefore at the edges 
of the crystals and in the vicinity of the boundary between two crystals the pit density 
is greater. А group of etch pits are enlarged in Fig. 13, and the border between the 
two crystals is shown enlarged in Fig. 14. In addition to the dislocation pits on the 
surfaces subjected to slow etching, fine figures arranged in rows of a density of 10%/cm2 
were observed. They may reflect the fine stresses (Figs. 14 and 15). 

In Fig. 16 the central part represents lines of a character similar to the “twisted 
cord“ observed on the cleavage plane (111). On the etched natural faces of the crystal 
of orientation (111) there also appear etch pits. Their size varies from 0.02 to 0.1 mm. 
The variation in the size of the large pits may be the result of their different density. 
At large dislocation densities sinall etch pits arise because of the cancelling of stresses 
resulting from single dislocations. The larger etch pits arise at the place of single 
dislocations and produce large local stresses (Read 1957). 

In some places on the etched natural faces (111) there appear pits shaped like - 
incompletely formed hexahedrons (Fig. 19). 

Etched surfaces close to the (100) plane were sparsely covered with pits of a 
completely different shape than on the (111) surface (Fig. 20). 

Besides the etch pits, melting patterns (observed by Millea and Tomizuka 1956) 
are visible (Figs. 20 and 21). The boundary composed of a number of pits lying along 
one line may be seen in Fig. 22. 

Pits situated very close to one another are frequently deformed (Figs. 23 and 24). 
It is possible that they reflect different dislocations and the mutual effect on one another 
causes, changes in the crystal lattice deformation and hence in the shape of the pits 
formed by etching. | 
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On surfaces making large angles with the planes (111) and (110) no etch pits are 


observed. The purity of the samples investigated corresponds to a conductance hole 
concentration of 1016 cm”? (p type). 


Conclusions 


l. The shape of etch pits depends on the crystallographic direction. Surfaces 
making small angles with the planes (111) are covered by pits similar to those observed 
on germanium crystals by other authors. This is most probably caused by the similar 
crystallographic structure. Surfaces making small angles with the planes (110) are cove- 
red with pits indicating uniform or greatly restricted local deformation of the lattice. 

2. In the case of slow etching of the (111) surfaces over a long period of time at 
depths which exclude any influence due to grinding and polishing, there were observed 
fine figures with a density of about 1000 times greater than the large etch pits occurring 
simultaneously. 

3. The large pits disclosed by etching of the planes (110) and (111) probably 
correspond to single dislocations; the fine figures may correspond to small periodic 
stresses resulting from the interaction between dislocations. 

4. The “twisted cords observed on the plane (111) may be either a system of 
single screw dislocations with large Burgers' vectors (giant dislocations) or small angle 


Fig. 25. (a) Fig. 26. 


boundary dislocations. A possible mechanism for the formation of dislocations in 
the first case is shown in Figs. 25 and 26. Fig. 26a represents the system of planes 
before the displacement, while Fig. 26b shows the same set after the slipping dis- 
placement. | ? 

The authors would like to express their indebtedness to Professor L. Sosnowski 
for his valuable remarks in the course of this work and to Mr. W. Giriat for supplying 
the indium antimonide crystals grown by him. The authors also thank Mr. J. Rautusz- 
kiewicz for the determinations of degree of crystal purity. 


List of Etching Agents Employed 


No. 1 (rapid etching) No. 2 (slow etching) 
HNO 50% HNO, 30% 
HF А 50% HF 20% 


CH,COOH 50% 
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No. 3 (slow etching) No. 4 (rapid etching) 
HNO, 40% HNO, 50% 
HF 3308 HF 50% 
CH¿OH 25% with small hydroquinone 
impurity 


КРАТКОЕ СОДЕРЖАНИЕ 


Я. Аулейтнер, Б. Колаковский. Микроскоповые наблюдения поверхностей моно- 
кристаллов InSb. 


В этой работе представлено результаты наблюдения некотрых поверхностей 
монокристаллов InSb травленных и вылущенных. Иеследовано поверхности 
с кристаллографическими индексами (110) и (111) и другие. Установлено, что на 
разных кристаллографических поверхностях форма ямок травления является 
разной. Наблюдано два типа фигур травления отличающихся друг от друга по 
величине и плотности. Для выявления фигур подобрано соответствующие тра- 
вящие средства. 
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RADIATION ACCOMPANYING PULSES IN 
A GEIGER-MULLER COUNTER 


By Kazimierz W. OsTROWSKI AND LESZEK TUREK 
Institute of General Physics of the Mining and Metallurgy Academy, Cracow 
(Received September 17, 1957) 


The ultraviolet radiation in a Geiger-Miiller counter was investigated by means of 
photo-sensitive counters with quartz windows facing one another and connected to a de- 
layed coincidence circuit. Two components of this radiation were distinguished, a stronger 
component occurring within the limits of 0—2 psec (prompt radiation) and at least one 
order of magnitude weaker than a second compo nent occurring with a delay of 2—300 usec 
(delayed radiation). 

The prompt radiation has a continuous spe ctrum with a maximum at about 2700 A. 
This radiation is interpreted as radiation from the recombination of argon ions. 

The delayed radiation is strongly depen dent on the counter supply voltage and its 
intensity can be measured only beginning from about 70 V before the end of the plateau. 
The maximum intensity of the delayed radiation when the counter is normally filled with 
an argon and alcohol mixture occurs at about 50 usec. The delay corresponding to the 
radiation maximum for a higher argon pressure decreases with an increase in pressure. 
The delayed radiation is interpreted in accordance with the Wiser and Krumbein hypo- 
thesis as originating from the metastable argon molecules. 


Introduction 


Investigations of spurious pulses in Geiger-Miiller counters thus far have been 

based on the analysis of the observed time intervals between pulses of the same counter 

` аз determined by means of delayed coincidences. The measurements are greatly 
` facilitated by exposing the counter to pulses of radiation at regular intervals of time. 
In experiments of this kind two maxima in the time distribution of the delayed coin- 
cidences have been observed (Montgomerry and Villard 1947, Putman 1948, Fenton 
and Fuller 1949, Curan and Rae 1947, Kupperian, Murray and Feeny 1951, 
Müller 1953, and Puppi and Gill 1955). The first maximum was interpreted as 
mpurities (electro-negative gases 0, H,O) of the 


originating in the negative ions of i : | | 
= Е filled. The second is associated with the process 


mixture with which the counter was 


of neutralization of the ions at the cathodes. | 
The above methods do not explain the source of the photons which produce 


| photoelectrons by means of a photo-effect at the cathode during the counter discharge. 
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If electro-negative gases are present, these photoelectrons may produce negative 
ions which cause spurious counts. The unsuitability of this method is based on the 
necessity of using “bad counters in the measurements, and therefore counters 
not working stably. Such counters are used to heighten the effects being mea- 
sured. This method thus gives a distorted picture of the operating mechanism of 
a good counter. In addition, the results of the above-mentioned investigations refer 
only to the problem of the slope of the characteristic of a “bad“ counter and do 
not explain the mechanism limiting the plateau length. 

Purely experimental work in connection with plateau length was performed by 
Weisz (1952) on the effect of the pressure of ether, by Rochester and McCusker (1945) 
on the effect of the pressure of argon and alcohol on the plateau length, and Barrére 
(1951) on the effect of the purity of xenon on the plateau length. The paper of Wiser 
and Krumbein (1955) put forward the basic idea that the formation of metastable 
atoms in the argon (*P,) state is closely connected with the ending of the plateau; 
in three-body collisions these may form metastable molecules of argon (A,*) of rather 
considerable stability. 

In case of spontaneous decay of the A, molecules, or in the case of their collision 
with alcohol molecules, radiation is emitted. This radiation may produce a spurious 
count if it is emitted after an interval of time longer than the dead time of the counter. 
A mechanism for the formation of A,* molecules from atoms in the metastable state 
(ЗР,) in three-body collisions was given by Coli (1954), according to whom the third 

2005 argon atom (150) is necessary in this phenomenon to 
conserve momentum and energy. The dependence of 
the mean time required for the formation of the argon 
molecule on the square of the argon pressure sup- 
ports this ‘suggestion. 

An idea similar to that of Wiser and Krumbein 
is found in the paper by Höhle (1956) for low-voltage 
counters using a helium and neon mixture. 

Particularly suitable for the experimental inves- 
tigation of ultraviolet radiation in G—M counters as 
discussed by Wiser and Krumbein, and Höhle are 

: “а t: the photosensitive counters with a platinum cathode 
Fig. 1. Principle of the radiation : : > 
ее ВО. used in the present work. Their photoelectric effi- 


20 cm 


1mg Ra Pb 
. 10cm 


f= source counter ciency and design is given in an article by Buja et 
D — detector counter al. (1957). 
Equipment 


Fig. 1 shows the principle on which the measurement of the G-M counter radiation 
is based. Counter S, hereafter called the source counter, is irradiated by a beam of 
gamma rays from 1 mg Ra. Counter S as well as detector counter D have 
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quartz windows facing one another. The detector counter should be as sensitive as 
possible to ultraviolet radiation. Counters with platinum electrodes were used here. 
The measurements were made by two coincidence units. Unit 1 had a resolving 
time of 1.3—2.0 usec and contained a delay line which gave an adjustable delay in 
steps of 0.7 usec within the limits of 0 to 7.7 usec. The diagram of the coincidence 
unit I is shown in Fig. 2. 

Unit II had a resolving time about 10 usec and permitted delays adjustable within 
the limits of 10 to 500 usec (the delay was given by a multivibrator with an adjustable 
pulse width) (Curran and Rae 1947). The coincidences were counted by means of 
a scalar. The counting rate of both counters were measuered separately by means of 
an ratemeter. 


JI) 


Fig. 3. Spectrometer for analyzing the radiation spectrum of G—M counters 


A quartz spectrometer with a high sensitivity was used to investigate the spectral 
distribution. The lens diameter was 3.5 cm, and the focal length of the collimator 
was 7.5 cm, while that of the objective was 25 cm. The source-counter wire was set in 
the place of the collimator slit. The receiving slit in front of the movable counter 
(detector) had a width of 2 mm. The counter (detector) position was controlled by 
a screw fitted with a micrometer scale. The detector counter and radioactive preparation 
were shielded by lead blocks (Fig. 3). An Original Hanau S 300 mercury lamp with 
a suitable attenuator was used to calibrate the spectrometer. 

For the sake of certainty filters (cellophane, glass 0.1 mm, lucite) were also used. 
Their transmission was determined on the spectrograph with the help of the photo- 
sensitive counter. 


Results 
Visual observation and detection by means of photographic plates 


With the eye accomodated to darkness and with a counting rate of the order 
of 10° cpm the glow about the wire of operating counter was plainly visible. For a vol- 
tage of about 300 V above the threshold even the single flashes corresponding to the 
. individual pulses could be seen. A distinct darkening of the photographic plate placed 
directly into the counter took place in the course of one minute at a counting rate 
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o 5 
of the order oF ех 10° cpm. Attempts to use a photographic plate to determine the 
spectral distribution of the radiation with the use of a Soviet ISP 22 spectrograph 
and even a spectrometer of greater sensitivity were not successful. 


Use of photo-sensitive counters to detect radiation emitted by a counter 


In a simple arrangement such as that shown in Fig. 1 the counter radiation may 
be readily recorded by a counter connected to a scaler or integrator. If, for example, 
the source counter 5 records 2 х 10% cpm and the background of the detector counter 
D due to the scattering of gamma rays is 2 x 103 cpm, then the measured effect is 
anywhere up to ten to twenty thousand cpm. 

Measurements of delayed coincidences between pulses from counters S and D 
indicated that the source S emits at least 90% of the photons in 0—2 usec after com- 
mencement of the discharge inside it. Closer investigations showed that the ultraviolet 
radiation emitted during the counter discharge can be separated into two components, 
one of greater intensity emitted in 0—2 usec after passage of the quantum producing 
the pulse in the counter — henceforth called the prompt radiation, and a component 
of smaller intensity emitted with a delay greater than 2 usec — henceforth called 
delayed radiation. 


Properties of the prompt radiation 


Fig. 4 is a plot of the intensity of the prompt radiation vs the overvoltage of 
counter S. Measurements with different cathode materials and counters of different 
diameter showed that this radiation, calculated for some solid angle, does not depend 
on the previously mentioned factors and that it is not absorbed significantly by the 
alcohol vapour. Investigation of the absorption of this radiation with different materials 
used as filters showed that this radiation is not monochromatic. 
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Fig. 4. Plot of the prompt radiation intensity vs the overvoltage. The counting rate cf the source counter 
. 4. S 
is 1.5 x 10% cpm. The detector counter counts up to 60% of the source counter rate’. 
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Fig. 5 shows the spectrum of the prompt radiation from a counter filled with 
an argon and alcohol mixture and measured in a spectromter of high light sensitivity 
specially built for this purpose. In Fig. 4 the difference in dispersion for different 
wavelengths was not taken into account. At the maximum of the spectrum 120 coin- 
cidences per minute were recorded. Counter 5 counted about 50.000 cpm. 
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Fig. 5. Prompt radiation spectrum. The ordinate axis represents a quantity proportional to the number 

of quanta emitted in the interval of wavelengths covered by the slit. The counting rate of the source counter 

was 50000 cpm. The voltage is 200 V above threshold. There were 120 coincid./min recorded at the 
maximum, 


Similar results were obtained for an argon-methylal counter, but the results in 
this case are less accurate because of the poorer stability of the small-diameter methylal 
counter and the high counting rates used in this experiment. 


Properties of the delayed radiation 


Small delays. Fig. 6 shows the time distribution of the delayed coincidences 
measured with coincidence unit I. It turns out that the intensity of the delayed radiation 
is strongly dependent on the overvoltage to a far greater degree than in the case of 
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Fig. 6. Time distribution of delayed coincidences for an overvoltage of 200 V. The counter was filled 
with 90 mm Hg of argon and 10 mm Hg of ethyl alcohol. The source counter pere 18 000 cpm, 
the detector counter 11 000 cpm. 
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prompt radiation. In a normal counter, i. e. a counter whose plateau lenth is about 
250 V, the radiation for small delays can be seen with an overvoltage greater than 
180 V. The intensity of the radiation at a small delay of 7.7 usec was 20 coincid./min 
with counter S registering 18 000 cpm and counter D 10 000 cpm, primarily due to 
prompt radiation. The number of accidental and Compton coincidences was of the 
orderof 1 coincid./min. 

Similar results were obtained for the counter filled with the argon and methylal 
mixture. | 

Similar results, even more distinct, were obtained when a single counter with 
the anode divided by beads was used. The counting rate employed was 10° cpm owing 
to which the effect was more marked. In this case the delayed coincidences for pulses 
from the same counter, but connected to the separate channels in the coincidence 


unit, was investigated. 


Large delay (above 10 usec) 


Figs 7 and 8 represent the coincidence rate for large delays. With an increase in 
the delay we observe the coincidence rate first to increase and then to fall. The position 
of the maximum coincidence rate depends on the pressure of the argon in the counter. 

Fig. 9 is a plot of the intensity of radiation for a delay of 150 usec as a function 


of the overvoltage (curve A). For comparison, a similar curve for the prompt radiation 


is given for the same conditions (curve B). Curves similar to curve A were obtained 
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Fig. 7. Time distribution of delayed coincidences for overvoltages of 280 V (B) ang 300 V (A). The 
RE was filled with 90 mm Hg ofargon and 10 mm Hg of ethyl alcohol. The counting rate of the source 
counter was 20 000 cpm. 


for a delay of 7.7 изес for counters with diameters of 14 and 25 mm filled with an 
argon and alcohol mixture or an argon and methylal mixture. This decreases the possi- 


bility of confusing the delayed radiation with the prompt ай sana by a pulse 
after the counter dead time. At a counting rate of 104 — 6 x 10% cpm no 
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maximum was observed as should be expected to be the case after the dead timè. This 
is because the effect of the delayed radiation under the conditions of measurement, 
in particular, near the end of the plateau, is more distinct. In order to find the maximum 
associated with the dead time it is necessary to operate the counter near the beginning 
of the plateau and increase the efficiency of the registration (window dimensions, 
efficiency of the detector counter). 
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j 
Fig. 8. Time distribution of delayed coincid ences for an overvoltage of 250 V. The argon pressure was 
190 mm Hg, the ethyl alcohol pressure was 10 mm Hg. The half-life for the delayed radiation decay is 
estimated as 30 usec on the basis of Figs. 7 and 8. As may be seen from Figs. 7 and 8, the delay for the 
radiation maximum is clearly dependent on the pressure. 


It did not prove possible to investigate the radiation spectrum. Investigations 
using filters indicated that the spectrum has a complex structure. One basic difficulty 
in method used thus far to record the delayed radiation should be noted. Near the 
end of the plateau the counting rates of both counters are close to one another. In 
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Fig. 9. Comparison of the dependence of intensity of delayed radiation (Curve A) and prompt radiation 
(Curve B) on voltage. $ 


such a situation the detector counter, because of its dead time, loses part of the pulses 
from the delayed radiation. Improvement of the geometry of the arrangement by 
increasing size of the counter window does not greatly improve the efficiency of 
registration. The authors avoided this difficulty by using counters with a long plateau. 
The plateau length can be relatively easily increased by coating the interior of the 
counter cathode with China ink. This effect is explained below. 
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Discussion 


The authors interpret the prompt radiation to be the limiting continuum (Grenz- 
kontinuum) of the recombination of argon ions to the levels of excited argon. The 
radiation from recombination to the ground state 15, corresponding to the very 
strongly absorbed vacuum ultraviolet cannot be taken into account in our consider- 
ations since the alcohol vapour and the quartz window certainly absorb it. If in 
accordance with the literature we take the recombination factor equal to about 2 x 10-19 
(Kenty 1928), we obtain a good value for the order of the recombination radiation. 
А more accurate comparison cannot be made because of the unknown ratio of the 
number of argon ions to that of alcohol in the first phase of the discharge, the neglect 
of the recombination to the ground state of argon, and the effect of the electron 
temperature on the recombination factor. 

The spectral distribution with a maximum of about 2700 À is related to the 
distribution of the kinetic energy of the electrons and the dependence of the probability 
for the recombination of ions on the energy of the electrons captured by them. The 
energy of the electrons in recombining is added to the energy difference between the 
ionization (15.7 eV) and excited levels. If it is assumed that we have to do primarily 
with the first four excited levels 1P,, 3P,, *P,, 3P, (11.5 — 11.8 eV) and that the 
probality of recombination does not depend too strongly on the electron energy, then 
the electron temperature may be estimated from the wavelength corresponding to 
the maximum of the prompt radiation spectrum. The expected long wave cut-off of 
the recombination radiation is 15.7 — 11.8 — 3.9 eV, i. e. about 3200 À, the mean 
electron energy is 4.6 — 4.0 — 0.6 eV, where 4.6 eV corresponds to the wavelength 
of 2700 À and 4.0 — 15.7 — 11.7. The electron temperature about the counter wire 
(anode) is thus about 7500°. 

The fact that the dependence of the prompt radiation intensity on the over- 
voltage is stronger than a linear one is explained by the electrons ionizing chiefly the 
alcohol molecules, which have a lower ionization potential. As the voltage increases 
the number of argon atoms formed increases in relation to the total number of ions. 

The authors put forward the following arguments against the suggestion that 
` the observed radiation is the radiation which, as Korff predicted long ago, accompanies 
the charge exchange between the argon ion and argon molecule: 

1. The spectrum of this radiation is not linear as should be expected in the case 
that the above prediction is correct (4.4 eV, i..e. about 2800 A). Quanta of energy Z 


4.4 eV are present. | $ Я 
2. The radiation in question is emitted in a time corresponding to the time for 


collecting the electrons by the anode. 

Any eventual doubts relating to the phenomenon of ion exchange can be settled 
by further measurements of the spectra, in particular, measurements of the dependence 
` on the voltage and investigation of the phenomenon in other quenching gases and 


non-self-quenching counters. . 
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The prompt radiation plays a basic role from the standpoint of the quality of the 
counter in two cases: 

1. When the counter contains electro-negative gases 

2. in a beaded counter. 

Ad 2. As a result of the prompt radiation in a beaded counter the discharge 
spreads to some extent from one section to the other by means of a photo-effect in 
its cathode; this leads to an increase in the effective resolving time of the counter. 
This phenomenon is more distinct the higher the voltage that is applied to the counter 
and the higher its photo-sensitivity. The photo-sensitivity of the counter increases 
with aging (see e. g. Jurkiewicz and Gierula 1951). As the resolving time changes 
with the aging of the counter the accuracy of the counter measurements at high count- 
ing rates decreases. In addition, in a beaded counter the spreading of the discharge 
from section to section by means of the delayed radiation should be taken into account, 
since this limits the plateau length. There still exists the phenomenon of the dependence 
of photo-sensitivity on the counting rate. This causes a further complication. 

It is because of these effects that these counters are not more widely used. 

The prompt radiation has already been measured independently of our work by 
Huster and Ziegler (1956). 

Delayed Radiation. This radiation is a confirmation of the predictions of Wiser 
and Krumbein (1955). According to these authors, in three-body collisions with argon 
atoms, the metastable atoms of argon (3P, 11.7, 3P, 11.5 eV), as already mentioned 
previously, give metastable argon molecules Az. Their spontaneous decay gives 
delayed radiation. 

Parallelly to the above-mentioned processes, the extinction of the metastable 
states of the argon atoms and the argon molecule takes place by collision with argon 
and alcohol particles. 

Wiser and Krumbein propose that these collisions accompany the radiation, 
where, as a result of the predissociation effect, radiation accompanies approx. every 
million decay collisions. Delayed radiation caused mainly by spontaneous decay of 
metastable argon molecules can produce spurious pulses if they are emitted after . 
a time exceeding the dead time and give a photo-electron. Basically, the dependence 
of the delayed radiation intensity on the overvoltage has a similar character as in the 
case of prompt radiation (Fig. 9.). We may account for slight difference between curves 
A and B in the following way: the probability for the recombination to metastable 
states 3P, and 3P, assumedly is strongly dependent on the electron temperature 
(counter operating voltage) as for the remaining states. In addition, part of the meta- 
stable states of the argon atoms arise directly from collisions with an electron and not 
through recombination. 

According to Coli (1954) the half-life for the decay of the argon molecule is 3.4 usec 
and does not depend on the pressure; the energy of the photons radiated during the 
decay is about 10 eV (1200 À). Our results indicate that this time is of an order of 
magnitude greater; we record photons directly in the region 2000—3000 À. For 


. 


Radiation Accompanying Pulses in a G-M Counter 107 


metastable argon atoms Coli gives a half-life of T = 1/9 p?, where p is the pressure 
in mm H. This formula does not stand in contradiction with our data. 

The discrepancy in the results as regards the half-life of the decay of the argon 
molecule may result e. g. from the existence of the metastable argon states 3P, (11.5 eV) 
and 3P, (11.7 eV). These states can give metastable molecules of different stability. 
Because of the use of argon as a scintillator gas (Forte 1956) a final explanation of 
the discrepancy in the results is of interest. 

It is known from theory and experiment that the counter dead time suddenly 
decreases above a certain voltage (about 200 V above threshold; the constant M 
derived by Wilkenson (1948) expressing the ratio of the produced charge to the charge 
associated with the wire capacity is then equal to 1). Thus, because of the exponential 
decrease of the delayed radiation in time, the mechanism leading to the ending of the 
plateau is then explained. 

Let us try to explain othér fasts by means of the above-described mechanism. The 
increase in argon pressure increases the counter dead-time. On the one hand, the 
number of metastable molecules formed should increase, since the probability of 
triple collisions will be proportional to the square of the argon pressure; on the other 
hand it should decrease, since the free path of the electron decreases, which favours 
the formation of alcohol ions over the excited states of argon. The net effect is that 
we observe an increase in plateau length as the argon pressure increases and then 
a drop (Rochester and McCusker 1945), which are the results of the present 
authors. 

Similarly, it is long known that the alcohol pressure favourably influences the 
plateau length up to fifteen mm Hg. This is connected with the extinction of the 
A, metastable states by alcohol. 

' The influence of the counter diameter on the plateau length is also known. This 
is connected with the dead time. In order to obtain good counters with small diameters 
it is helpful to increase the argon pressure. | 

Finally, the diminishing of the photo-sensitivity (e. g. by coating with carbon 
or by using a cathode of the Maze type) has a marked effect on the plateau length. 
In this way counters may be obtained with a plateau length of about 400 V. This effect 
was employed on a large scale in the production of glass counters with a graphitized 
cathode and in counters of the Maze type. By increasing the counter: diameter, e. g. 
to 50 mm, a large increase in plateau length cannot be obtained. This is related to the 
dependence of the dead time on the voltage. At about 250 V above threshold the 
dead time begins to decrease rapidly with the voltage. наш, the plateau length of 
about 250 V is, in a way, the “normal plateau length*. ca 

A plot of ‘the delayed radiation intensity vs the voltage perans tne учи. на: 
of the change in the counter plateau length to be determined if its active lengt © 
hanged by some definite factor. For example, let us assume that 


hoto-sensitivity is с - 
bis produced counters have a plateau of 250 V. What will be the length of plateau 


for counters that are twice as long? To find out we construct a curve with ordinates 
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twice the size of the ordinates of the curve representing the variation of the delayed 
radiation with voltage (Fig. 10). 

We next produce a horizontal line from point A to the curve 2 l and we obtain 
a new plateau length 20 V smaller. Since the changes in dead time are not taken into 
account, this procedure gives only the upper limit for changes when the counter is 
lengthened, and the lower limit when the counter is shortened. 
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Fig. 10. Graphic method for calculating the limits of the changes in plateau length produced by changes 
in the active length of the counter or changes in photo-sensitivity. The ordinate axis is in arbitrary units. 


At high counting rates (105 — 108 cpm) the resolving time decreases even to 
0.25 times the value measured at low counting rates. A good counter, nevertheless, 
is stable under these conditions if the voltage is not too near the end of the plateau. 

This is in agreement with our results, namely the delayed radiation normally 
begins at about 170 V above threshold. 


Conclusions 


It seems that the method presented for investigating the phenomena occurring 
during the discharge of G—M counters by investigating the radiation arising during 
different phases of the discharge is most appropriate if one is interested in the details 
of the mechanism of counter operation. It seems that the method described can be 
` used for investigating certain special effects such as the recombination of ions and 
metastable states of the gas used for filling the counter. Perhaps after some improve- 
ments this method can be used to investigate in detail the phenomenon of neutraliza- 
tion of ions at the cathode. 

According to the review article by Massey (1955), these phenomena have not yet 
been sufficiently investigated. Also, they are important for both G—M and scintillation 
gas counters. 
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КРАТКОЕ СОДЕРЖАНИЕ 


К. Островский и Л. Турек. Ультрафиолетовое излучение в счётчике Г.М. 


Пользуясь фоточувствительными счётчиками исследовано ультрафиолетовое 
изучение счётчика Г. М. Счётчики были приложены к себе кварцевыми окош- 
ками и соединенные с цепью опозданных совпадении. Обнаружено, что излуче- 
ние состоит из двух компонент: первой в пределе от 0—2 сек. (незамедленное 
излучение) и второй на один порядок величины меньшей, которая появляется 
с замедлением о 2—300 писек. (замедленное излучение). : 

Спектр незамедленного излучения непрерывный, с максимумом около 2700 А. 
Излучение объясняется, как излучение рекомбинации йонов аргона. 

Замедленное излучение очень строго зависит от напряжения на счётчике 
и только начиная с 70 V от конца плато можно его регистрировать. Максимум 
замедленного излучения для нормального наполнения смесью аргона и алкоголя 
появляется при замедлении около 50 исек. Для больших давлении аргона заме- 
дление, которое отвечает максимуму падает с увеличением давления. 

Происхождение замедленного излучения от мэтастабильных молекул аргона, 
объясняется согласно гипотезе Висэра и Крумбеина. 
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DETERMINATION OF THE ORIENTATION OF MOSAIC BLOCKS 
BY MEANS OF A FINE-FOCUS X-RAY TUBE 


By JULIAN AULEYTNER 
Institute of Physics, Polish Academy of Sciences and Institute of Physics, Warsaw University 


( Received September 26, 1957) 


In the method described in this paper, an X-ray tube with a point focus, constructed 
by the author, was used as a source of X-rays. From the point focus with an effective 
diameter of the order of 40u the radiation falls on a narrow slit between two cylinders. 
In this way the chopped-off beam falls on the surface of the crystal under investigation. 
Two photographs are made of a suitable diffraction line Cug,: (1) with the crystal 
oscillating and the photographic film immobile and (2) with the motion of the film coupled 
to the motion of the crystal. 

In this way the disorientation of the blocks of the mosaic can be determined for 
angles greater than 30°” and, in particular cases, of twist of the order of seconds. The method 
permits the obtaining of a picture of the mosaic structure over a larger area of monocrystal 
surface in a relatively short time and the disclosure of the direct dependence of the Bragg 
line intensity on the mosaic structure. 


In the majority of cases real crystals have a mosaic structure, that is, they consist 
of blocks of microscopic size and twisted with respect to one another by very small 
angles (of the order of seconds or minutes). 

To investigate the disorientation of the blocks of the mosaic, X-ray spectrometry 
with a single, double, and even triple crystal and ionization methods of recording the 
diffracted beams are used (James 1950, page 280; Renninger 1955; Backovsky 1956). 
Photographic recording is rarely used because of necessity of employing monochrom- 
atization of the X-ray beam by means of monocrystals and because the moving of the 
photographic film during irradiation requires relatively long exposures (Intrater et 
al. 1954, Weissman et al. 1951, Weissman 1956). 

In the method described below, an X-ray tube with a point focus constructed 
by the author (Auleytner 1957) has been used as the source of X-rays. The diffracted 
beam is recorded by means of a photographic film. F ig. 1 shows a diagram of the 
apparatus and the path of the radiation. From the point focus S with an effective 
diameter of the order of 40 y. the X-ray beam falls on a narrow slit between two cylinders. 
The diameter of the cylinders is 8 mm. The distances between the focus and the slit and 
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between the slit and the crystal are 60 and 40 mm, respectively. The beam chopped off 
in this manner falls on the surface of the crystal under investigation. The path of the 
beam coincides with the axis of rotation of the crystal. The crystal may oscillate about 
the position satisfying the Bragg condition for the interplanar distance associated with 
the surface being investigated and for a given wavelength of the characteristic spectrum. 
The diffraction pattern is recorded on an immobile plate F located at a distance of 
100 mm from place where the initial beam falls. In the case of a crystal with a block 


Fig. 1. 


structure the Bragg diffraction line is not of a homogeneous compositions. As a result 
of the twisting of the blocks as well as local damage to the surface there appears either 
a reinforcing or a weakening along the lines (Figs. 2a and 3a). If during the irradiation 
the photographic film oscillates with the crystal, then the individual blocks satisfying 
the Bragg condition for the same crystallographic direction at different times give 
the diffraction pattern recorded at different points of the plate. 

In Fig. 2 the Bragg lines obtained for the cleavage plane (100) of a NaCl crystal 
can be compared for two cases: (a) an immobile film (b) the film oscillating with the 
crystal. In Fig. 2b may be seen the displacement of the individual elements from which 
the line shown in Fig. 2a is composed. The individual elements of line 2a indicate 
a twisting of the macroblocks within the limits of from 0.5’ to 15’. The gaps visible 
on line 2b correspond to a strongly twisted fine block. й 

Fig. За and 3b compare the diffraction pattern obtained by Ше same method for 
the case of KBr monocrystals subjected fo the action of nitrogen oxides. Fig. 3b indi- 
cates the pressure of a continuous change of twist of the fine blocks of the mosaic 
within the limits of from 1’ to 25’ (in the plane perpendicular to the axis of rotation 
of the crystal) and a strong deformation of the entire crystal. A discontinuous darken- 
ing corresponding to the individual blocks of the mosaic is visible in the direction of 
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the diffraction line. Figs. 4a and 4b show the diffraction pattern obtained in the same 
way for a quartz monocrystal slab. The slab surface was matted. A microphotogram 
of the K, Cu line was obtained for the case in which the crystal and photographic 
film oscillated in unison, giving rise to a line that is diffuse in comparison to that 
obtained when only the crystal oscillated (Fig. 5a). The diffuse line may be caused by 
the superposition of effects resulting from a small disorientation of submicroscopic 
blocks of the mosaic and a divergence of the beam falling on the crystal. The photo- 
graphing of the beam falling right at the slit and at the place where the crystal is 
mounted permitted an accurate determination of the beam divergence. At a slit width 
of 40 y the beam divergence was 12”. Thus the diffuse line observed in case 5b is the 
result of the mosaic structure of the crystal. 

The use of a slit width of 10 u permits the obtaining of a diffraction line 30” wide. 

Fig. 9 shows a microscope photograph (magnification 73 X ) of the Ka, and 
Ka, lines of copper obtained by coupling the oscillations of the crystal to the oscillations 
of the photographic film in the case of Ge (11) plane, dislocation density 10°/cm?, 
smallest angular width of the line 34”). 

By shifting the crystal and photographing each time the diffracted beam on the 
oscillating film, we can obtain a picture of the structure of the entire surface in which 
we are interested. The exposure time for a slit 40 p wide was always 15 min. for 
oscillations of the film and crystal of the order of 10°. By limiting this range to 2° the 
exposure time could be shortened 5-fold. In the experiments, the position of the table 
should be fixed in such a way that the axis of the rotation of the crystal lies precisely 


on the surface under investigation and coincides with the plane of the'incident beam. 


A special sighting tube was used for this purpose.. 

The diffraction lines are so narrow that photometry does not always accurately 
reflect the shape. For this reason microscope photographs at a magnification of 75 X 
should also be used for line analysis. For example the photograph in Fig. 6b obtained 
for a NaCl crystal can be used. At the points marked by arrows the line of Fig. 6b 
is so sharp that the lines corresponding to Ka, Cu and Ka, Cu are distinguishable. 
This part of the line is shown in a microscope photograph at a magnification of 73 times 
(Fig. 7). The smallest angular width of the Ka, line standing out on the diffuse back- 
ground is 1. The diffuse region has an angular size greater than Г, and, in example 
cited, it results from the twist -of the fine blocks of the mosaic. 

A microradiogram of the NaCl crystal, obtained by irradiating the surface with 
a wide beam of radiation coming from the point focus and by photographing the 
Bragg diffraction line confirms the existence of a block structure in the NaCl crystal 
(Fig. 8). The points marked by arrows correspond to the deformed parts of the line 
Е prc of the described method is testified to by the fact that it permits the 
direct disclosure of the influence of the crystal's mosaic structure on the intensity of 


the Bragg line. According to the theory of Darwin (confirmed by experiment, James 


1950) the intensity of the beam diffracted on a crystal surface of considerable thickness 


$ 
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is greater for a mosaic crystal than for an "ideal" one. Photographs of diffraction 
lines obtained for an optically smooth cleavage surface (100) of an NaCl crystal show 
different diffuseness in accordance with the degree of disorientation of the blocks 
(Fig. 6b). Points of greater darkening on the Bragg line (Fig. 6a) correspond to the 
greater diffuseness of the lines obtained by coupling the motion of the crystal to the 
motion of the photographic film (Fig. 6b). 

Further investigations concerning the limits of the possibilities and applications 
of this method will be made. 

Conclusions: (1) The method presented above permits the determination of the 
disorientation of blocks of the mosaic for angles greater than 30”; with such an arrange- 
ment of the blocks that the diffraction lines are displaced so that they lie one over 
the other (e. g. Fig. 2b) twists of the order of seconds can be determined. (2) By means 
of the method presented here a picture of the mosaic structure of a larger area of 
the monocrystal surface may be obtained in a relatively short time. (3) The method 
permits the disclosure of the dependence of the Bragg line intensity on the mosaic 
structure of the crystal. 

I would like to express my sincere gratitude to Professor L. Sosnowski for dis- 
cussion of the results of this work and for his constant interest in its progress. 


КРАТКОЕ СОДЕРЖАНИЕ 


Ю. Аулейтнер. Определение ориентации блоков мозаики с помощью острофокусной 
рентгеновской трубки. 


В описанном методе этой работы источником рентгеновских лучей являлась 
рентгеновская трубка с точечным фокусом, сделанная автором. Из точечного 
фокуса с эффективным сечением порядка 40 и, лучи падают на узкую щель 
между двумя цилиндрами. Вырезанный этим способом пучёк падает на поверх- 
ность исследуемого кристалла. Делаются два снимка соответствующей дифрак- 
ционной линии (нп. Сика): 1. когда кристалл осцилирует, а фотопластинка 
является неподвижной; 2. когда движение фотопластинки сопряжено с движе- 
нием кристалла. Этим способом можно определить дезориентировку блоков мо- 
заики для углов больше 30”, а в исключительных случаях скручивания порядка 
секунд. Метод позвает получить в относительно короткое время образ мозаиковой 
структуры больших областей поверхности монокристалла, а также непосред- 


ственное обнаружение зависимости интенсивности линии Брагга от мозаико- 
вого строения. | 
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ON THE MODIFICATION OF THE FEYNMAN’S 
"INTEGRAL — OVER — ALL — PATHS* METHOD 
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Feynman's “sum-over-all-paths“ method was already successfully applied for 
solving problems of the fundamental equations of quantum mechanics, of the theory of 
stochastic processes, as well as of sum-over-states. In this note we give a modification 
of the mathematical formalism of the Feynman's method as well as an approximate 
solution of the problems mentioned above for large "time" intervals and a discrete 
spectrum of states. 


1. Introduction 


We begin with the rigorous solution of some equations of the diffusion type. 
Further, we investigate Feynman's “integral - over — all — paths“ and use our 
previous results to modify the computational technique of this method. One may 
say roughly that the relation between our methód of summation and that used by 
Feynman (1948) is similar to the relation between the Lebensgue'ian and Riemannian 
integrals. The outcome of some considerations of a combinatorial and probabilistic 
character is a formula approximating with great accuracy the solution of the diffusion 
problem mentioned above for large values of the t-parameter. 

A shorter note on the same subject will appear in Il Nuovo Cimento. 


_2. Rigorous solution of equations of the diffusion type 


2.1. Let us call an equation of the form 


aphid _ 
ot 


= lim RE (j |j) di' v Gs t 49 — Y (js; t—4 1 (1) 


41>0 
The unknown function is Y (7, Jo» 2), and z(i|k) 


an equation of the diffusion type. | 20 
s n matrix for a t-interval A ¢; the initial 


are given functions; they form a transitio 


f (117) 
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condition for £ — 0 is 


pee a pce me 194 Ја) ШЛ is continuous, 
APS ra if j is discrete. 2) 


JJo 


To this type of equations belong e. g. the Schródinger equation, equations of stochastic 
processes, as well as equations describing the sum-over-states (in all these examples, 
except for the last, t has the meaning of time). If the j-spectrum of states is not continu- 
ous we replace the integral on the right-hand side of (1) by a sum. 

For the sake of simplicity we confine ourselves in our considerations to the case 
in which the j-spectrum is discrete and the number of levels is bounded, say 1. We 
assume also that the variable t is discrete. Finally, we will discuss the case of continuous 
t. Thus we consider the equation 


l 
Фут) = 2 2 GI) YO’ оп — 1), (la) 


SL 


Where lL т апа j, Na Wl s n а 
2. 2. To solve (la) we define the generating function for y (7, Jo; п): 


Ф (j, jo; 1) = 2 a" y (js jo; n) (3) 


uniformly convergent in the neighbourhood of A = 01. It follows from (la) that Ф 
satisfies the equation 


Ф (j, jo) = ûj, +4 21261») Ф (у... (4) 
y 
Equation (4) yields the solution 
aa ode 
Ф (12702 A) > y, 9 
where 
_| —Az(1|1) +1,—Az(2] )),... — åz (1 |1) 
za [2201222012 +1.  —4=(12 
—Az(1|l),—Az(2|J),... —Az(l|l) +1 


is the determinant of the “stochastic matrix (we write here inverted commas because 
z (i | Е) must not be probabilities) and J,, is the minor of J (A). We put this solution 
into the form 


2029 = sO = (5a) 


1 * = v 
If we need to use the continuous t-parameter, we must replace the power series on the right- 


hand side of (3) by a Laplace transform of y (0). 
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with 


0)) = — E In J (A). (5b) 


It is obvious from (3) and (5a) that, if we develop Q in a power series in respect to A 


Q (4) = za c (п). (6a) 
we get 
ee до (п 
y join) = CI. (6b) 


The solution (5) or (6) has an interesting feature: each solution of a problem 
concerning a process of diffusion (Markoff processes, Schródinger equation, etc.) 
утау be written as a partial derivative with respect to a transition matrix element, 
say 2 (хо | x), of a“potential function“ о (t). This function is symmetrical in z (7 | k). 
For example, we are able to write all elements of the S—Matrix as “gradient compon- 
ents of only one such “potential function" in the z (i | k)-space. 
2.3. We now prove the relation 


ER eio Aa HII «6| =p (1) 
J з 2 ¡=1 k=1 IH II н! 
і Е 


where the summation is subject to the following restrictions: 


1 
Sy ee (7a) 


Let us write the right-hand side of (7) in the form 


n=0 пъ 4] С, С] 
1 d ; | Lae Lo 
«ПУ: | П я: 6) 
k=1 Li=1 k=1 


where we sum over м» for which 


X щеп +1 (8а) 
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holds, and the integrals are taken along closed paths C; in Ше ¢;— complex-planes 
enclosing the origins of the coordinates. We are now able to perform in (8) the sum- 
mation over n. We have 


I 92 ЫР z (i 954 Ars = 
(ва) hu па 0) 


We introduce new variables 


than we have 


consequently 


(a) ф- jn 15095 5 
з fis (10) 


where the paths C'; are similar to C;, but in the 2, — planes. 
2.4. The function 


JESA KA (11) 


tends to zero together with A. Consequently 
20) = тва + ак) 


may be developed in the еч of A = 0 in a power series in A. In this way 
we get 


n+1 


a, 


s=1 j=1 


7. 
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where 


1 1 x 
ПП +—— =Y°... a, 12 
2 k=1i=1 П I ux! 2 2 Ins «e Вр» (12а) 


the summation on the left-hand side is subject to condition (7a). 


3. On Feynman's “integral — over — all — paths“ 


3. 1. Feynman (1948) proposed a new method for solving the Schródinger 
equation by means of so called "integral-over-all-paths". 

In this section we will present a modification of this method; we do not, 
however, ‘have in mind the modification of the main idea of this method but only 
the mathematical formalism. 

Let us illustrate (i) the method of Feynman and (i i) our approach by the example 
of the Schródinger equation in one-dimensional space?. 

We seek the solution 


вол) = < аер I ие» 


of the equation 


200 
Y PEEL урен (13) 


2 
Ot 2m da? 


th 
with initial condition 
y (ль! = 0) = д(х— у. (14) 


This solution may be written by means of Feynman’s method 


y (x, Xo; t) = Ша (а sl dx, dX, 5425 


„0 —— m 
an п-1 


¡EA EU aj. 2 15a 
х I (кў) e» (xis 4-1) (15a) 
where 
| | 1 Рта | 
y e E At ) — И (xi) ) (15b) 


| 3 in a one-dimensional space; the generalization to more than one 


2 We perform our computation 
dimension does not seem to present difficulty and is outlined in section 4.2 of this note. 


у 


\ 


122 J. Lopuszanski 


and x, = x. In (15a) 


(16) 


тан exp E L (m)! = € exp I-i нац ха» 


i.e. Ше transition probability amplitude from state x; to state x; in Ше time 
interval At. If we denote, as in Section 2, the transition amplitude from state Е to j 
in the time interval At by z (j |k), the solution (15) may, of course, be written 


y (5, 29; f) = lim |... ан H z (x; | n 1). (150) 


At ^0 س‎ 
nát—t п-1 


One concludes from (15c) and Section 2 that Feynman's method may be applied 
to each equation of the diffusion type; the solution of such an equation has the form 
(15c) and may be immediately written down if only we know the transition matrix 
z (j| Е). 

So far Feynman’s idea. 

An (n — 1)—fold integration appears in formulae (15а); this summation pro- 
cedure was already modified by some authors (see e. g. Davison 1954). We propose, 
next, another way of computation: we replace the (n — 1)—fold integration by sums 
over groups of “similar“ paths. What will be understood by “similar paths“ will 
be clarified further. on. 

3. 2. Let us demonstrate our attempt by the simplest example, namely, when x 
has discrete levels and the number of levels is bounded, as in Section 2. Although 
this case is meaningless for physical practice (except perhaps the case of spin waves 
in a crystal), nevertheless, this example is suitable for the clear presentation of the 
main idea of our speculations, without troublesome computations. This case was 
already investigated and rigorously solved in Section 2. As in Section 2, we also assume 
that time is a discrete variable, the shortest interval being A t. 


Thus the solution (15c) becomes 
1 


1 п 
У. У d sx | (15d) 


ж=1  xy-1=1 ¡=1 


This solution consists of a sum over “chains* of length п, which all start in x, and 
end in x. Such a “chain may be also expressed by | 


Й «6s ICI EAM | (17) 


j=1 k=1 


where m,, denotes how many times the term г ( | Е) appears in the “chain* concerned. 


We have, of course, 
7 


чет 
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n 
безл. е. с , d : Е : Ч Е = 
Each “chain“ [Į z (x;|x; ,) may be written in a unique manner in the form (17) 


{= 
on the right-hand side, but not vice versa; the number of “chains“ which may be 


characterized by the same ensemble of exponents m;, and which begin in x, and end 
in x we denote by 


<x | M o (19) 


we quote by М the. dependence of (19) on all matrix elements ту. We mention 
that if M is given, the length п may immediately be fixed by means of (18a). 
The necessary condition to form at least one “chain“, if М is fixed, is 


> ma, = > пу for kA x, xo (18b) 
i PE 

and 2 m, + 1 -D SUD (18c) 
j j 
Ут, —1 = ту. (18d) 
7 7 


Relations (18b—d) simply mean that in a “chain* Ше number of transitions towards 
a state, say k, must be equal to the number of transitions away from this state, x and хо 
excepted. If we exclude from our considerations the case in which M is only diagonal, 
then (18) is also a sufficient condition to build a “chain“ of length п. It follows at 
once from (18) that 


У ZUA m= a toi m 


consequently the final state x is uniquely determined by M and хо. Thus we call some 
“chains“ or paths similar" if they are characterized by the same М and xp. 
Using notation (17) and (19), the solution (15d) may be written 


(id NA 
w(x, xg; n) = Dy П П=0[" <x | M | x>, (15e) 
М j=1k=1 
the summation is extended only to such М?» for which (18a) holds. Of course, these 
<x | M | хо? °з, for which (18b—d) do not hold, vanish. | 
The expression оп the right-hand side of (15е) has the form of a multiple power 
series of the variable г (j | Е); thus the (x | М |.x9)’s play the réle of power coefficients. 
It should be emphasized that the nature of the physical problem is entirely in- 
volved in the 2( | Е) (е. с. the mode of interaction of the particles of the system is 
described by their potential energy appearing only in the exponent of the паар 
probability amplitude, similarly to the example of the one-dimensional Schródinger 
equation in the x-representation, already discussed in this section). On the other 
hand, the expressions <x | М | хо» do not at all depend on the physical nature of the 
problem; to find them we must only employ geometrical and combinatorial consider- 


ations. 
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The problem of evaluation of the <x | М | хо в may be called “the problem of 
fence construction” (see fig. 1), Indeed, let us imagine we have a transport of planks 
of A t cm width. Each of the planks is painted blue on one side and red on the other. 
From among these planks m;, are cut off slantwise so that, looking at the red side, 
the left border has the height k cm and the right —7 cm (see fig. 2). The total number 


ee <. 
At 
----4..- 
Fig. 1. Fig. 2. 


of planks is n. Imagine, further, that we have two poles at a distance n At cm one 
from another. One has the height x, cm, the other — x cm. 

We ask in how many ways one may build a fence between these poles, so that 
the upper edges of the planks form a continuous, broken line, the ends of which are 
of the same height as the poles (see fig. 1). Each side of the fence is uniform in colour; 
if we look at the red side we should have the pole of height x, cm on our left. 

We now show how we can evaluate approximately, though with great accuracy, 
the number of fences, i. e. the expression (x | M | x,» for large п. To do so we use 
combinatorics and make extensive use of the law of large numbers. j 

З. З. Let us take u, planks of width At cm representing transitions (k — j) 
(БЕ = 1, 2,...1), satisfying (7а). We try to construct from these planks a fence in the 
form of a circle (see fig. 3), the upper sides of the planks form a continuous, broken 
line. We are, of course, able to construct at least one such a fence because of 


> 216 — B ua = 0, 


provided that џ,, 0 for 15 k. à 

Assume that we know all possible modes of constructing such a round fence. 
Among these modes some of them are "essentially different“, i. e. that no two fences 
become identical by rotation. The number of such essentially different fences we 
denote by C (и), where и denotes the matrix Hine 

Let us now consider one particular round fence. From this fence we are able, 
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in general, to construct several fences of length п At cm, shown in fig. 1. For this 
purpose we remove from this round fence a plank of type (x => xp). Since in one 
round fence the plank of type (х > хо) appears Uxx times, consequently we can obtain 
from one round fence at most My» different fences of length n At cm. In the over- 
whelming majority of cases we get exactly и. fences. The exceptions come from 
those round fences of which some parts appear repeatedly; thus it may happen that 


by removing a plank of type (x — xp) in two different places we get the same fence 
of length n At cm. 


If we neglect these exceptional cases, we can write 


(x | M | хо? = Uxx C (A), (21) 
where, in accordance with (7a) and (18b—d), 
Mir zi бу, д» + Ше С 
Substituting (21) in (15e), we obtain 


д 
Е [1 [120| me ceo). (22) 


the summation being restricted by (7a). 
3. 4. Taking into account the fact that С (и) does not depend on x and хо and 
is symmetrical in all ру, we may surely make use of the law of large numbers for 
the purpose of determining С (и), without risking the removal of the dependence 
of x and x. 
The law of large numbers reads 
— — P (Mp (il) 


lim P | = e) = 1, (23a) 


п со 


where Р (k) denotes the probability of finding the system in the state k, p(i | А) is the 
conditional probability of finding it in the state i, if it was, before in the state k, and 
в > 0. In other words, the “chains” characterized by mj, which differ slightly from 


ma = n P (k) p (i | B) (23b) 


are very probable. The remaining chains taken together give a very small contribution 
to the probability. We note that the number of “highly probable chains“ is usually 
very small in comparison to the number of “not-very-probable chains“. 

3 To show this, let us Т the special probabilistic scheme р (i| ®Ю = 1/1, (i, k = 1,2, ... D); then 
each chain has the probability =” and 


o n 
mik = Г 


According to the law of large numbers, 
P» TI Пр" <2z|M|%>=1 
мм: R 
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From (23) it follows that the function 
П Пр (i | ky" <x | M | xo >, (24) 
i k 


considered as a function in m;,-space, almost always has an exceedingly pronounced 
maximum in the region given by (23b), 1. e. 
0 
Mik 


=P (i | А). (23c) 


> Mir 
7 


We must take into account the fact that this maximum may be found from (24) under 
the supplementary conditions (7a). We notice also that our considerations are suitable · 


for large п, consequently we may employ и» with the same effect, instead of mj, 
as well as 


П П p (i | K)"i* C (и), (25) 


instead of (24). 
The usual Lagrangian multiplier method yields the result 


д x ; 
a, mew) Амр В +4 4- 4j — Ar = 0. (26) 
Lyk 
This equation should yield us the values m;,, provided that the function С (и) is known. 
This, however, is not the case. Although we do not know С (м), we know approximately 
the values of the maximum, i. e. (23c). Substituting from (23c) the value for p (i | k) 


or 


l-n Y) <2|M|m>=1 
M~M 


Since the total number of chains is J”, thus > <х|М | хо > exhausts almost all chains; here, the 


MM 
dependence on x, for such a rough calculation is meaningless. It follows that, if we take another probabi- 
listic scheme, where 


y 1 
р (i| Б) = те? 


Ше maximum of 


П Пос| та «а | Мо» 


i k 


. . n 
will be, of course, different from ту, = 4 „then the chains characterized by m;, lying near this maximum 


form a very poor class with respect to the class of chains characterized by mj, ~ я 5 
Thus, the case p(i | А) = 1/1 is the only exception when the class of “highly probable chains‘ is 
very numerous. : 


+ de 
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in (26) we obtain differential equations for С(и), puc when solved, yields 


У} иж a 
ОРТ М 
k=? Пе 


year 


(27) 


These, of course, are merely approximate equations, proper only for such и which 


may be a maximum for a certain probabilistic scheme р (i | k); e. g. they must satisfy 
the relation 

Min = Щр. 
Nevertheless, we are able to extrapolate (27) for all м satisfying (7a). In (27) 


A = const : exp [— A (n + 1)}, 


because 
exp {— > 2 Ay) и} = 1 
j 
in virtue of (7а). Taking into account the fact that u,, are large numbers, we may write 


T 
C(u) — A {ка (28) 
k=1 Пи! 


To evaluate A, we use 
1 
У Xs, бш) =, (29) 
у=1 p 


where the summation is taken over и, satisfying (Та). With help of (28), we rewrite 
(29) in the form 


SE UI fof Engr 


Ml 
j=l Lik * 


RATE 
PIA Yu) $ fon p NU "LU 
where Uu; = Mig — ды дк and consequently ш; = W— дъ. We have, further, 


ATE Lele) 3 f Ye) I = 
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iFx 
ON ur 
пада П 1 и 1)! (uy — 1)! 
т=Ех,у 
Mr)! n! = | 
+An М NA RC е ... Ус A (n+1) Ir, 
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Thus 
== | (п +1), (30) 


and, finally, with help of (22), (28), and (30), we obtain 


5 1 П! 
EIA A у н A я 
_ O9z(xy|x) [n+1 m ПП 15 ^5 ПП» (31) 


Y (Ху Хо 5 n) 


where Ше summation is restricted by (Та). One sees that the expression in the brackets 
in (31) is the first term in the alternating series of w (n) in (12) divided by the number 
of terms in this series. 

Looking at (7) we find immediately that 
5 со 


Ф (x, pd a 


OU) = 


"Nes 
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We now check if the function satisfies equation (4) approximately. Divergences 
appear only for small n. 


Since D (x, Xo pA) = 
fg (ткт 1 $- f’ dtr д dl, 
Жар 271 TD A + 6651» Glo Ge 
we have 
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1 
(x | y) =) 4 dep d$, сх 
X е 
XL (i |») E A =) }- {У ll Se 


x E Е ) | de | (33) 
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The second term in the brackets of (23) may be written 


utum $ $0 бйз 1 > 
[E x; сша 4210 &—54@|%)% = 
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because only the terms belonging to large n are of importance, and thus we may replace 


IS 


by 


СЭ Eorp 


The first term in the brackets of (33) vanishes fof х Æ ху, because 
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since 
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1 
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1 1 hood of С, = 0. 
is a one-valued function of ¿, in the neighbour 
Thus, by (34) and (35), we see from (33) that equation (4) is satisfied for x A Хо: 
For x = x, the first term in the brackets of (23) yields 


E 1 
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Thus the right-hand side of (33) may be written, by (7), 


= n+1 
i за [гу => a 


za Шия! 
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It is evident that for large п and sufficiently large u,, the last term may be neglected 
in comparison to the first. Consequently equation (4) is approximately satisfied for 
X = ха too. 


4. Transition to the continuous t—variable 


4. 1. We have shown that the exact solution of an equation of the diffusion type 
is given by (12) and an approximate solution, proper for large n, has the form (31). 
All this was done for the case in which the # — variable takes discrete values. Let 
us now investigate the case in which A tends to zero. 

We define the function p(s), s — 0,1; p(0) — b; p(1) — a; a and b are arbi- 
trary complex numbers. Further, we introduce by a linear transformation a new 
variable 


s(a +b)—a 

Vab 1 
We write p(s) — q(£). It may be shown, as usual, with the help. of characteristic 
functions (fe Appendix), that for 


р, = У... 210 (Ê) 9 (5), | (37) 


Е = (36) 


where the summation is taken -over &, satisfying 


Е = a (E, + te Hs £,), (37a) 
Va 
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a central limit theorem is valid; indeed, for large n, 


(a E b)r+1 (a + p)? n : 2 
T скри ey ҳ ; 
OE EE: V27 nab P| 2n ab (X a): ws 


= 


in other words, 
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We apply this result to the computations in (31) for At— 0. We confine our- 
selves to the case in which the number of levels is bounded, or, at least, the spectrum 
concerned is discrete; the investigations of ihe case of continuous spectrum meet 
some obstacles. 


The most suitable case for ореол is that of stochastic processes, because 
all the z (j | k) are > 0 and (a + b) = х: (j | k) = 1. If we consider Ше sum-over- 


-states some complications already arise, ia ра condition z (¡| Е) > 0 holds 
too. The most troublesome case is, however, that of quantum mechanics, because 
z (j| Е) are complex numbers‘, Let us pay some attention to this case. 

By (39), for large u, we are able to replace 


———— ПИ 
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4 Since 
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and 


and since the solution with respect to the #-- variable is continuous, we have фот each representation 


1 
У): (0 =1+ © (A 1), 
j=1 


Hence, as in case of probabilities, the normalization condition holds for At — 0. 
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where Q-! is the reciprocal quadratic form to 
О (%1... 21) = > »4 Е Yo Ont? y YY X ar is (41a) 
r,=0 r,=0 r= 1=0 


and o? is the determinant of the matrix of "dispersion" 9,, y, ++» "1-1 belonging to 
the “distribution function“ p(s,, 59, ... 5j 4) 3 S1, За: 5-1 = 0,1, (see Appendix): 


p(0, 0, ...0) = z (1| #) (X z (j |)? =b, 
р (1,0, ...0) = z (114) (1:010) =a 


p (0,1, 0) = z (2 |А) (12010) = ae (41b) 


p (0,0, ...1) =z (1—1 |4) (2010) = ai; 


otherwise = 0. 
The essential feature of (41) is that it has the form of a manifold Gaussian distribution 
and that the dispersion is proportional to и, in the exponent. If we go over to the 
continuous t—variable we must carry out the transformation 


ljk At = typ; fra, 118-0 
м At tg. (42) 


Thus we obtain a new Gaussian distribution, but with a dispersion equal to t, o? At, 
proportional to A t. It is obvious that, if A t tends to zero, these distributions will 
tend to a Dirac delta function. We shall not, however, replace these expressions directly 
by Dirac delta’s; we see that the main contribution in (31) comes from the region? 


(43) 


Since >) дь = p мк» it follows immediately from (43) that these t, must satisfy the 
7 7 ; 


equations 
| | 
t; съ Y tp zll А | (44) 


which are already the equations for the stationary state 


t 


p (j; > со) ©; 


5 This becomes plausible if, for rough computations, we employ for (41) the method of steepest 
decent in the case of real functions, or the method of stationary phase in the case of complex functions. 
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from these equations we can find tj. Thus, of interest are the contributions due to 
groups of paths, which show a deviation with respect to the stochastically typical 
group of paths, the latter leading to stationary state. Computations of these deviations 
need some care and will be the subject of a separate paper; the starting points for 


these computations will be equations (12) and (31) as well as (41). Thus in the future 
we shall have to deal with 


j=1 


1 Ke 5 = 
n +1 Уп(У 019) Ола) 3 сар 
Т 


Е E 1 d tr z (1 | Е) 
> т, y (45) 


where the summation over Т, i. e. over t; is restricted by 


4. 2. We mention, finally, that the idea presented in this note may be immediately 
generalized from one to several dimensions. For this purpose the scalar indices ] and 
k in 2(j | k) as well as in m,, must be replaced by vector indices, e. g.. 


j Ts Uv Jo: Jp k = (hy, ko, ә. kp); 


consequently we have also 


z (j | К) and my. 


It is, however, obvious that the numerical computations will be performed most 
easily in the case of one dimension. И 

If, in addition, only the elements 2 (j |j) and z (j |J + 1) are different from zero, 
which is approximately realized in several physically important cases, the matrix ид 
reads, е. g. for хо < x 


Bi T Ма for b= ло and pde 


ДАЙГА = Ша +} for xy <k < a. 


This means that it is almost a symmetrical matrix; consequently the principle of 
detailed balance is here satisfied. In this approximation the solution yields only the 
stationary part. Since и; = >) Hie the number of independent variables reduces 
* 1 
k 


to [. 
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Appendix 


We shall show that for P,(&) defined by (37) the central limit theorem (38) is 
true. 


To show this, we take the Fourier transform of both sides of (37). We obtain 


t n 
= (0) = |" (а) | 
үл 
where we use the notation 


Pa (1) = У exp {it E) P, (E); 


the summation is taken over 


(A 
¡=1 b 
where s; = 0,1, and 
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t== u ры tE == 
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+ 
Let us develop q, (t) in a power series in t, which is always possible; thus 


Pn (t) = (a + Б)" ( o se SI 


where the terms of higher order involve powers of n^! higher than the first. For n tend- 
ing to infinity we get 


: Pn (t) -5. 
in (297) ver 


Now, the inverse Fourier transform yields 


1 53 
Р, (E) ~ (a +b)" dé * 


=e 
2л 
Since ё is given by (А) and consequently 
ise d +6 
Vnab 
then (38) follows immediately. 


In Ше many-dimensional case, i. e. for p(s}, Sa, ... 51-1), 5; = 0, 1, defined by 
(41b), we proceed similarly as in the case of | — 2; we have here 
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where OO бра) em p (81... 511) and the summation is taken over ES, which 
satisfies 


РЕ 
4=1 


related to s? by 


a A а; 
Ei = s — == Ss ОЩЕ 
2 a; +b 
7 
We carry out the multi-dimensional Fourier transformation with respect to 
INL: 
ty И-» с 
qat. tre 91 | Ee ES Е 
EET. 
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For n-> со we obtain the asymptotic formula 
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The inverse Fourier transform of g, reads YA, = 7 
j 
1—1 


pan a? Yate er uw (Vn &, ... Vn &-1) а 


j=1 


where the notation used here has been given earlier. We notice that Q^! is pro- 


portional to o. 


КРАТКОЕ СОДЕРЖАНИЕ 


Я. Лопушанский. Модификация метода Фейнмэна „интеграла по всем траек- 
ториям“. 

Мы даем строгое решение некоторых уравнений диффузионного типа. Далее 
мы рассматриваем Фейнмэновый интервал „по всем траекториям” и излагаем 
новый математический подход к методу Фейнмэна; для этой цели мы исполь- 
зовали полученное раньше нами решение. Упращая можно сказать, что соотно- 
шение между нашим и Фейнмэновым методом произведения суммирования по- 
хоже соотношению между интегралом по Лебеду и по Риману. Мы даем дальше 
приближенное решение уравнений типа диффузии, применимое для больших 
промежутков времени и дискретного спектра состояний. Вид этого решения 
относительно несложный. 
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Conservation of the third component of iso-momentum (angular momentum in 
the space of isotopic spin) can be secured even for weak interactions at the cost of assuming 
a composite structure of particles in isospace (two independent isospins J and K). Particles 
may be classified according to quantum numbers b, K?, Ks, I, І, where b is the baryon 
number. This classification enables to arrange particles into a table exhibiting interesting 
regularities. Rest masses are explicable in terms of the above quantum numbers. There 
exists a fine structure and a hyperfine structure of the mass levels described by J? and J. 

Anomalous behaviour of the spin K (two orientations only) suggests isospace to be 
a euclidean subspace of a ps-euclidean space. It is suggested to combine isospace with 
ordinary space into a more-dimensional ps-euclidean space. Then the anomalous behaviour 
of the spin K may be interpreted in terms of a Zitterbewegung. 

Assuming a 6 or 12-dimensional space it is easy to understand why pion is 


pseudoscalar. 


Introduction 


Let us begin with a few questions: 

1. What is hidden behind the group-theoretical classification of interactions: 
(i) strong interactions characterized by a full symmetry with respect to the group 
of rotations in three dimensions,.(ii) electromagnetic interactions having an axial 
symmetry (symmetry with respect to rotations in two dimensions) and (iii) weak 
interactions spoiling even an axial symmetry? In particular: 

la. Is it possible to say that there exists an isospace similarly as one uses to say 

at ordinary space exists? 

= lb. If * ще how to explain the privileged role of a “third” axis? 

le. Are weak interactions really incompatible with the conservation law for the 


- third component of the angular iso-momentum (isospin) ? | | 
ld. Is it to be regarded for established, that leptons have nothing to do with 


isospin and isospace ? 
(137) 
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2. Is it possible to derive, from purely group-theoretical considerations, some 
information about the masses of elementary particles without taking advantage of 
the computations of self energies 


3. Is there any connection between isospace and ordinary space? 


4. Why Nature seems to prefer field quantities with certain transformation 
properties as compared with others, apparently equally well possible? In particular: 
why is pion ps-scalar and not scalar? 


We shall try in this paper to give a common answer to all these questions which, 
at the first sight, seem to be disconnected or only loosely connected with each other. 


The great advantages brought about by the conception of isospin in the problem 
of nuclear forces, in particular the laws of conservation of 1? and Да, strongly suggest 
that isospin is not a purely formal coriception but that there is an isospace with a si- 
milar amount of reality as the ordinary space. On the other hand, the violation of 
the law of conservation of 1? by electromagnetic interactions and of 1, by weak inter- 
actions suggest just the opposite. 


If there is a violation of symmetry in the ordinary space, we are able to explain 
this fact by some external circumstances (e. g. by the existence of an external magnetic 
field). Such circumstances are principally removable. But the violations of symmetry 
in isospace do not seem to depend on any external and removable circumstances 
but seem to exhibit an intrinsic anisotropy of isospace itself. A conception of an 
anisotropic space, without pointing out any excuse for it, seems inadmissible. Thus, 
stressing the anisotropic features brought about by electromagnetic and weak inter- 
actions, we are rather inclined to regard isospace as fictitious. 


The above dilemma could be solved in favour of a physical reality of isospace if 
we succeeded in demonstrating that the violation of symmetry is only apparent, or 
if we succeeded in understanding the reasons of the deviation from a full symmetry. 


First of all we shall show that the violation of the third component of iso-momen- 
tum by weak interactions may be regarded as apparent only. This will reestablish the 
axial symmetry and constitute a new argument in favour of the existence of an isospace. 
The question why there is a privileged third axis will be discussed later, together with 
the problem of relation between isospace and ordinary space. 


1. Conservation of the “third component“ of angular isomomentum 


As is well known there exists a relation between the third component of isospin 
I; and the electric charge. The simplest way of writting this relation is q +I, + x = 
= 0 where x = 0 in the case of the pion, the hyperon A, the hyperon Х, and for 
their antiparticles, whereas x = 1/2 for the heavy meson K** the nucleon N, and 


the hyperon Е, but x = — 1/2 for K~*, the antinucleon N and the antihyperon =. 
The quantity x is additive. Strong as well as electromagnetic interactions conserve 
the total I; of the system, so that also the total x is conserved. For weak interactions 
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among the above mentioned particles (to be called henceforth hot particles) there is 
a change A I, = + 1/2. Since q is conserved, this means at the same time a change 
Ax = + 12. Thus, x plays the same role as strangeness and the selection rule of 
Gell-Mann (1954) and Nishijima (1954) may be stated as follows: 

For strong and electromagnetic interactions A x = 0, for weak interactions 
A x = + 1/2 whereas no processes occur with dx = +1. 


It is strongly suggestive to regard x as a third component of another isospin J and 
to write 


—q=1,+J; (1) 


This formula implies that the electric charge q is a result of two kinds of rotation in 
isospace around a third axis, and that the third component of the total angular iso- 
momentum is always conserved. 

The appearance of two distinct iso-momenta means that elementary particles are 
in general composite structures in isospace, somewhat similar to the atom (in the 


ordinary space). The particle consists of an “iso-nucleus“ with spin J and of a "shell" 


with isospin E The total angular iso-momentum of the particle consists of these two 
parts. Its third component is always conserved but the difference between strong and 
electromagnetic interactions on one hand and weak interactions on the other hand 
consists in that the former couple only nuclei with nuclei and shells with shells whereas 
the latter couple nuclei with shells so that there is an exchange of momentum A I; = 

= — ДЛ, = +} among them. i 

In this way the law of conservation of the third component of angular isospin 
has been generalized for the case of weak interactions among hot particles. The problem 
of extending these conceptions for leptons will be dealt with later. 

With the aid of (1) the hot particles can be arranged into the table 1. According 
to this table К+» together with K — may be regarded as a single particle existing in 
four states. Also N together with & may be regarded as a single particle. Since I and 
J are simultaneously integer or half integer, the particles possess always an integer 
total angular isospin and may be called “iso-bosons”. 


140 J. Rayski 


2. A classification of elementary particles 


The characterization of particles by means of two types of quantum numbers 
I2, Iş, J?, J4 has the advantage of a great simplicity but has an unsatisfactory feature: 
pion is characterized by exactly the same quantum numbers as У, and Х by the same 
as Из antiparticle. Also the К meson has exactly the same quantum numbers as the 
system JV, Я (as well as their antiparticles). The hyperon A has the same quantum 
numbers as its antiparticle. Therefore the above characterization is obviously insuffic- 
ient. The lack of uniqueness of characterization is due to the fact that we have 
disregarded a very important quantum number b denoting the number of baryons: 
b — 1 for a baryon (nucleon or hyperon), b — —1 for an antibaryon and 5 — 0 
for a meson. Let us take account of b by generalizing the formula (1) 


b—q=1,+ К, (2) 


where K, is connected with J, by Кз = Ja + b. Owing to the law of conservation of 
the number of baryons, Кз is conserved for strong and electromagnetic interactions 
but is changed by ДК. = — AT, = + \ for weak interactions. Thus, Ка plays 
also the same role as strangeness. In contradistinction to Ла the new quantum number 
K, posesses an interesting property: it is positive for each particle and negative for 
each antiparticle. Arranging hot particles according to the increasing values of b, 
then according to the increasing values of Kg, then according to the increasing values 
of I, we get the table 2. (2; has been placed to the right of A since it is characterized by 
a larger value of [® than A). 


Table 2. 


We see that the new characterization introduces a one to one ‘correspondence be- 
tween different sets of quantum numbers and different types of particles. 

The table 2 exhibits interesting regularities resembling the periodic table of 
Mendelejew. In spite of the fact that we have arranged everything in a monotonic 
sequence of the numbers b, Кз, Is from the left to the right, there appears automati- 
cally a reflection symmetry with respect to the centre. The neutral pion is situated in 
the centre whereas the particles situated symmetrically to the right and left from the 


me) e Mo d 


` Let us assume that photon is sp 
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centre are charge conjugate one to the other (zx? is the only hot particle charge con- 


jugate to itself). All particles situated to the right from zx? are (or in the doubtful 
cases л”, K^"^, may be called) particles, those to the left are antiparticles. But the 
most astonishing fact is that the particles arranged themselves in an ascending sequence 
of masses according to their distance from the centre. This curious property holds 
not only for the big mass differences between the multiplets but also for the small 
mass differences within each charge multiplet: the charged pion is heavier than the 
neutral one, the charged K-meson is (probably) a little heavier than the neutral one, 
neutron is heavier than proton, hyperon 2” is heavier than 2” which, in turn, is 
heavier thanZ*. It may be expected that also Æ- will be proved heavier than Ж° in 
the future. Of course the particles situated symmetrically with respect to the origin 
possess equal masses. 

The regularities exhibited in the table 2 can hardly be regarded as purely acciden- 
tal. Most probably they are manifestations of some profound regularities in the particle 
structure. This structure was completely disregarded by the traditional quantum field 
theory. A future satisfactory theory of elementary particles will have to explain the 
regularities exhibited by the table 2 and to clarify the meaning of the formula (2). 

In contradistinction to the quantum number Ла introduced by (1) and assuming 
only the values 0, + Y, the new quantum number Ка assumes also the values + 1 
and + 3/,. Since b and q are strictly conserved, an interpretation of Kg as a third 
component of an isospin (angular iso-momentum) is desirable in order to preserve 
the conservation (of the third component) of the total angular iso-momentum Ig + 


+ Kg. However, we encounter here a difficulty: if the hyperons 4 and 2 possess an 


isospin К = /K/ = 1 then we should expect three values of Кз, i. e. we should expect 


the existence of the counterparts of the hyperons A and У with Ка = 0. Such particles 
do not exist. Therefore we have to swallow the fact that K is able only of two (instead 
of 2K +1) orientations: parallel and antiparallel to the privileged third axis. Such 
a restriction looks at the first sight very strange but, as will be shown in section 3, 
it implies important consequences regarding the nature of isospace. 

Assuming that K is an isospin that may be oriented only parallel and antiparallel 
to the third axis we see that N, N, E, E do not form a quadruplet with K — 3/2 but 


two different ‘doublets: N, N with K = 1/2, and Е, Я with К-+32, 
The assumption of only two orientations of one of the isospins Кз = + K enables 


to extend the conception of isospin also for the case of leptons. Let us assume in 


analogy with (2) 


where 7 and Г have to play similar roles for leptons as I and К did for hot particles, 


1. e. represent two kinds of isospins, the latter capable at most of two orientations. 
; inless (in isospace) but electron together with neutrino 


form a triplet i = 1, k = 0, whereas muon is a doublet i= 0, Е = 1. This yields 


# 
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a table 3 somewhat analogous to the table 2. Photon and neutrino are situated at the 
centre. The particles situated symmetrically with respect to the centre are charge 
conjugate one to the other and the masses increase with the distance from the centre. 
The law of conservation of the third component of the total angular isomomentum 
for the system composed of hot particles and leptons holds in the form J, + К. + 


+ 1з + kg = 6 — q. 


Table 3 


3. Ап interpretation of rest mass!) 


The formulae (2) and (3) and the Tables 2 and 3 show that the electric charge is 
connected with the third component and the mass with the absolute value of the angular 
momentum in isospace. 

The mass spectrum resembles somewhat the energy spectrum of the atom. The 
absolute value of b plays a role somewhat similar to the principal quantum number, 


and the quantity К = [К] = /K,/ describes a fine structure of the mass levels. Since 
b — 0 for the mesons but b = 1 for the baryons, the former are generally lighter 
than the latter. These two levels split according to the fine structure quantum number 
K. The mass of the heavy meson is larger than that of the pion because its quantum 
number K is larger. Similarly, the masses of baryons increase with K. Moreover, 
there are hyperfine structure effects connected with the “nuclear“ isospin Г. The 
main hyperfine effect is connected with РР, the small effect with Ig. The sign of the 
small mass correction is positive (negative) if the resultant K + Tis larger (smaller). 
In this way the question (2) in the Introduction has been answered affirmatively. 

Neglecting the hyperfine structure effects the mass levels can be described by 
the formula 


mi = Pl aq (4) 


TM e A O 
1 A little different interpretation has been published elsewhere (Rayski, 1957). 
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where / is a fundamental length (units с = h = 1 are used) and x(K) is a quantity 


playing the role of the length of the angular momentum К. In a classical limit x is 
identical with K but in quantum mechanics 


RE VA(K + 1) fon TA sees 
K+} for K=4, 8, ... 


Thus, for integer K the quantity x(K) is the length of the angular momentum vector 
but for half integer K it is the eigenvalue of the well known operator j of Dirac. This 
operator corresponds as well to the length of the angular momentum vector. The 
formula (4) together with (5) yields a vanishing mass for the pion which means that 
its whole mass is a hyperfine structure effect. The formula (4) takes also into account 
that the mass of the nucleon is about twice as large as that of the heavy meson. The 
ratio of the masses of thy hyperons A and X to the nucleonic mass is 2% which agrees 
exceedingly well with the experimental value in the case of A but agrees very poorly 
for X. However, the difference т; — m, may be regarded as a hyperfine structure 
effect (Т = 1 for X but I = 0 for A). Finally, the ratio of the mass of the hyperon 
= to the nucleonic mass is 2% which agrees very well with the corresponding experi- 
mental value. 


х( 


(5) 


The numerical agreement of the ratios of the A and Я masses to the 
nucleonic mass with the experimental values is so good that the fact of proportionality 
of the baryon masses to the root үх can be hardly considered as accidental, 

It is tempting to look for an “anschaulich* interpretation of this proportionality. 
The question is why baryon masses are proportional to the roots үх? Assume that 
the masses are determined by the spin K (whereas the baryonic number serves only 
to determine the characteristic baryonic length Г = mm = 2) Classically the 
spin K may be visualised with the help of a rotator model. Its angular momentum 
(in isospace) is classically х = pd where d is the rotator length and p is the momentum 
of the rotating point. The simplest and most symmetric assumption connecting the 
quantities p and d with one another and with the mass of the particle is m =p m d 
in natural units c = ^ = / = 1. This yields immediately m = Vx. The assumption 
p == d means that the rotator possesses a reciprocal symmetry. The assumption 
m — p (or m — d) exhibits also a symmetry between the usual Браво Nr mass 
plays its dynamical role) and the isospace (where the rotator 15 шло. 


zy ПЕНЯ between the radius of the rotator (representing the particle) and the particle 
mass resembles the situation in general relativity where the particle is represented bv a Schwarz- 
schield sphere with a radius proportional to the mass. This suggests that in А. пише sd n 
the isospace will be riemanian and particles will be represented by Schwarzschicid yia S oth, 
ordinary space and isospace,. Thus, particles possess internal degrees of RÎ in ка o ies: 
dinates on those spheres and their internal structure will be described by spheri я пет Б 
a description by means of spherical harmonics оп a sphere with radius proportional to the mass 1 


equivalent to our rotator model. 


$ 
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In spite of the fact that a real understanding of the connection of isospace with 
ordinary space, and therefore of the connection of isospace with mass, is still lacking, 
it is very satisfactory that the rest mass ceased to be a mere parameter but appears 
to be a function of some momenta. Rest mass is rest energy. It appears now similar 
to the kinetic energy since both are functions of some momenta. 

The appearance of the quantities of the type K(K +1) ог К + 1/2 in (4) and 
(5) constitutes an additional argument in favour of regarding K as an angular momen- 
tum. This argument, obtained from a qualitative and quantitative discussion of mass 
levels, is quite independent from the former (formerly we assumed K; to be an angular 
momentum component in order to save the law of conservation of the third component 
of angular momentum in isospace). Thus, both arguments strengthen each other 
and contribute to the opinion that isospace is a physical reality manifesting itself not 


only through the traditional isospin J and charge independence of nuclear forces 
but also throughout other properties of elementary particles, i. e. their rest mass and 


their spin K. 


4. A pseudo-euclidean character of isospace 


In the foregoing sections we encountered a difficulty that the angular momentum 
К possesses at most two orientations instead of 2K + 1. The question is whether it 
is possible to account for this restriction in a natural way? A similar situation is familiar 
within the usual space: the spin of the photon possesses only two orientations parallel 
and antiparallel to a privileged direction: its line of flight. The same is true for any 
particle traveling with a velocity of light. This is connected unavoidably with the 
ps-euclidean character of Minkowski’s space. Therefore the peculiar property of the 


spin K indicates that isospace must be of a similar nature as the ordinary space, namely 
it must be a three-dimensional, euclidean subspace of a more-dimensional ps-euclidean 
space. Thus, there exists (at least) one more dimension playing a similar role with 
respect to isospace as the time does with respect to the ordinary space. Let us call it 
“iso-time“. By introducing *iso-time* we get a possibility of motions in isospace 
with an “iso-velocity“. This constituent of the particle which carries the isomomentum 
K (the shell of the particle) moves in the direction of the third axis with the “iso-velo- 


city“ of light which explains why K must be oriented only parallel or antiparallel 
to the privileged axis. In this way the apparent difficulty with the interpretation of 


K as an isospin disappears. Instead, we > got some information about the structure of 
isospace. 

The above results are by no means unavoidably connected with the assumption 
of a quantity K, assuming larger values than №. The same conclusion can be reached 
also by limiting ourselves to the formula (1) (where Л, = 0, 14) provided nucleon and 
hyperon Е are not treated as a single particle (doublet with respect to Ла) but as two 
distinct particles. We find then that the spin J; of the nucleon can be only J; = — 15 
whereas that of the antinucleon can be only J = % (for the hyperon E the situation 
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is reversed). This situation is analogous to the neutrino of Lee and Yang (1957) in 
ordinary space. Nucleon appears to be left-handed and antinucleon right-handed 
in isospace. Such a situation is unavoidably connected with the ps-euclidean character 
of the underlying manifold. 

Assuming the particle model offered by the formula (2) the situation changes 
in so far that all particles acquire a right-handed зриаШу, all antiparticles a left- 


-handed spirality as regards the orientation of K with respect to the third axis. 


5. On a possible connection between isospace and Minkowski space 


Asis well known, there exists a possibility of regarding the positon as the negaton 
moving backward in time (Feynman 1949). However, by reversing the direction of 
time it is the negaton that appears as the positon moving backward in time. This 
shows that it is possible to connect the problem of charge conjugation with time 
reversal. A similar situation holds not only in the case of electron but also in the case 
of pion. An exchange of the roles of negative and positive pions can be achieved by 
means of time reversal and an exchange of the roles of bra and ket vectors whereby 
creation operators become annihilation operators and vice versa. But the same can 
be achieved by means of a rotation in isospace around the second axis through the 
angle л reversing the direction of the third axis. This shows a possible connection 
between the time reversal and the reversal of the third axis in isospace. It is suggestive 
to try an identification, 1. е. to look whether it is possible to regard the privileged 
third axis in isospace to be nothing else but the time axis. 

If such an identification is feasible, it will mean a fusion of a two-dimensional 
space of electric charge (t4 tə) with Minkowski space - (€, х1, xg, Хз) into a six-dimen- 
sional space (f, to, t3» X1. X2» хз) where t = із. This space is ps-euclidean with a signature 
+ + +>>—-—. Its time-like dimensions constitute a euclidean subspace (ty, ta, їз) 
to be interpreted as isospace. Of course, it must be checked carefuly whether such 
a fusion contains no contradiction with well established experimental facts. This will 
be examined in the next sections. For the present moment let us point out only that 
the pretended fusion seems promising: 

1. By amalgamating isospace with ordinary space into a more-dimensional ps- 
-euclidean space we fulfil the requirement of the foregoing section: isospace appears 
actually to be a euclidean subspace of a ps-euclidean space. This was necessary to 
account for the peculiar behaviour of the isospin K. — i 

2. The idea of fusing isospace with ordinary space into a six-dimensional space 
is attractive since it enables to incorporate isospin on equal footing with the ordinary 
spin. The group of rotations in our six-dimensional manifold повне two sue, 
ing subgroups of rotations in the euclidean subspaces (5; t2» ty) and Сащ ха). The 
operators of infinitesimal rotations in the latter define the mechanical angular momenta 
whereas those in the former have to be interpreted as the angular isomomenta, As 
` these two sorts of operators commute, we can introduce representations in a product 


4 
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form, i. e. attach to the field quantities two kinds of indices, one connected with iso- 
spins, the other with ordinary spins. 

3. This idea seems satisfactory from a broad viewpoint since it means a progress 
in the same direction as relativity. The theory of relativity extended the geometrical 
realm of physical phenomena by fusing apparently disconnected spaces: the ordinary 
three-dimensional space with a one-dimensional time. The same is endavoured now: 
two apparently disconnected spaces are fused into a single space whereby the third 
axis in isospace is amalgamated (identified) with the time axis. 

4. This idea enables to understand intuitively why antiparticles can be regarded 


as particles moving backward in time and why one of the isospins, namely K, can be 
oriented only parallel or antiparallel to a preferred axis. 


6. A twelve-dimensional space 


In the foregoing section it was suggested to amalgamate isospace with ordinary 
space into a ps-euclidean six-dimensional space. The question is not only whether 
such an extension is possible but also whether it is sufficient. We should not be affraid 
of the eventuality of increasing the number of dimensions as there are several hints 
that the manifold underlying physical phenomena is actually more-dimensional. 

First of all, there is an argument derived from’ the general theory of relativity: 
the riemannian space is curved and can be placed only in a more-dimensional flat 
space. The embedding space (with a necessary minimum of dimensions) should 
possess a physical meaning. It is usually believed that the embedding space is ten- 
-dimensional. 

A hint that the physical manifold is probably twelve-dimensional is furnished by 
the hyperon 2. It is a particle involving the largest number of states: three states of 
charze to be doubled by two states of spin again doubled because there are particles 
and antiparticles. In order to describe the hyperon Х by a single equation of a Dirac 
type we have to use y matrices with twelve rows and columns which suggests a twelve- 
-dimensional space. 

One of the possibilities is to regard the presumable twelve-dimensional space 
as a three-dimensional manifold of Minkowski spaces. If it is so, the total space possesses 
three real and nine imaginary axes (or vice versa). The time-like axes constitute a 
three-dimensional isospace. It is a subspace of a ps-euclidean space. 

In the realm of twelve dimensions one of Minkowski’s worlds is privileged. Its 
privileged role is to be regarded not ds an intrinsic anisotropy of the general twelve- 
dimensional space but as the anisotropy of the object placed in it: the totality of matter. 
Similarly as a given material particle with a non-vanishing rest mass distinguishes 
a system of reference in the Minkowski space, (namely its system of rest) which by 
no means invalidates the principle of relativity, similarly as a given photon distin- 
guishes a direction in the ordinary space (its direction of flight) whereby the isotropy 
of the ordinary space is not invalidated, also the totality of matter in the world dis- 
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tinguishes a certain Minkowski space, our proper world, whereby the objective sym- 
metry of the twelve-dimensional framework does not need to be invalidated. 

Let us denote the dimensions as follows: t, v, v vs, ty Wy Wa Шар t3 X1 Хо Ха where 
the last four coordinates mean our proper world. In the same sense as one of Minkowski 
spaces, also one of the axes in isospace (t, ty tg) is distinguished, namely, tz, and 
plays the role of our time. It is possible to develop a tensor calculus in twelve dimensions 
but, in view of the distinction introduced by the actual totality of matter, it is more 
convenient to limit ourselves to the description of physical phenomena in a system 
of reference that takes account of the actual situation (similarly as in a situation brought 
about by the existence of a magnetic field it is convenient to place a third axis in the 
direction of the field. Instead of twelve-tensors it will be more convenient to use 


a Characterization according to commuting subgroups of rotations in the twelve- 
-dimensional space. 


7. Why pion is ps-scalar? 


Let us consider the group of rotations in isospace. A set of field quantities trans- 
forming like a vector under rotations in (t, & #3) denotes a particle with isospin 1. How- 
ever, the field quantities describing the pion cannot be a vector in isospace since 
such a set has not a covariant meaning in the (proper) Minkowski’s world. This compli- 
cation is an immediate consequence of the identity of t4 with the time variable in our 
proper world. In order to have a field component фз transforming like a third compo- 
nent of a vector under rotations in (& ta tg) but like a scalar under rotations (in ts 2, 
хә хз) we have to assume that фз transforms like the product dt, dx, dx, ахз. For 
symmetry reasons let us assume that the three field quantities фу, Pa, Pg transform 
like products indicated in the table 4. This shows that Фа cannot be scalar in the 
proper Minkowski's world but must be pseudoscalar. 


Table 4 
(The case of isospin Г = 1) 
field quantity | 
transforms like dt, dv, dv, dv; | dtdw,dwdws dt; dx, dx; dx; 


The remaining components are ps-scalar in the other four-dimensional subspaces 
instead of in the proper Minkowski world. There is no contradiction in it since 


the signs of фу and pa, and therefore also their scalar or ps-scalar character, are im- 


material. Indeed, by performing a gauge transformation 
1 
prep where p= nm tim) (6) 


e reversed so that we can regard them as scalar or ps-scalar 


9 f d ar в 
шуде ч danger of contradiction. The only thing 


(in our proper world) as we like without any 
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that matters is the behaviour of фз under inversion of t = tg and of %1, ху, Хз, and this 
was shown to be unavoidably ps-scalar. We conclude that from the point of view of 
the Lorentz group pion cannot be scalar but must be ps-scalar. 

By reversing our argument we can say that scalar particles cannot form a triplet 
with the charge states +g, 0, —q. 

The possibility of demonstrating, purely theoretically, that pion must be ps- 
-scalar seems to constitute a crucial test of the correctness of the idea of unifying 
isospace (t; & tz) with the ordinary space (x, x, хз) by identifying the privileged third 
axis with the time. 


8. Why there are two kinds of isospin? 


The type of transformation indicated in the table 4 applies to particles possessing 
an isospin J = 1. However it is also possible to supply a field quantity with an index 
i denoting components of a true vector in isospace provided the orientation of the 
vector in the direction of the third axis is prevented (table 5). This happens if there 
is a movement in isospace with an iso-velocity of light in the direction of the privileged 
third axis. In this case we have a polarization perpendicular to the direction of flight. 


Table 5 
(The case of isospin K = 1) 


field quantity 


transforms like 


In the case of isospin K = 1 the vector of polarization possesses only two non-vanishing 
components in the directions of t, and fg. 

Question arises what is the meaning of a “motion in isospace with an isovelocity 
of light“? A general motion in our twelve-dimensional space with a velocity of light 
is characterized by 


3 3 

p dt? - ҳа v? + dw} + ал) (7) 
A motion in isospace means a reversal of the roles of time-like and space-like coordi- 
nates whereby the space-like coordinates play the role of “times and the time-like 
coordinates mean “position“, In particular, the motion in the direction of the third 
axis is characterized by 


3 
а = У, (dv? + dw? + ахд - (7) 


i=1 


which means a change of position in the course of the ordinary time. In particular, 


"e d 
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it can be a motion in the ordinary space 


dt? = 2j dx?. (24) 


It is seen that a motion in isospace in the direction of the third axis means (in particular) 
a motion in the ordinary space. Therefore we ask whether anything is known about 
the existence of such a motion having nothing to do with motion of the particle as 
a whole? Indeed, such a motion is known as Schródinger's Zitterbewegung. Therefore 
we conjecture that Zitterbewegung characterized by (7"), or a generalized Zitterbe- 
wegung characterized by (7’), are responsible for the peculiar behaviour of the isospin 
K. Conversly, since there is a possibility of a Zitterbewegung, the particle can possess 
two different kinds of isospin: besides the normal isospin I the particle can possess 
also ап *abnormal* isospin K. The normal isospin appears if the field quantities 
transform as indicated in the table 4, the abnormal appears if the field quantities 
transform as indicated in the table 5. Of course a particle can be also described by 
a tensor involving both kinds of iso-vector indices and can possess both kinds of 
isospin simultaneously. 

The above considerations can be extended for the case of half-integer isospins. 


9. Charge conjugation and. inversions 


From the point of view of a more-dimensional space time reversal is ambiguous. 
We can reverse time either by means of an inversion, or by means of a rotation (around 
the axis 1, through the angle л) in isospace. In the first case the handedness of the 
system of reference is changed, in the second case it 1s not. 

Let us define the operation of charge conjugation as an inversion of the second 


axis in іѕоѕрасе 
Ge |. (8) 
Charge conjugation (8) causes a change of sign of the third components of isospins. 


In order to handle easily the operation C we have to use a convenient representation. 
In the case of the pion it is convenient to use the representation 


( 1 А 
p= ФЕ, g! Bip agit uh P3- (9) 


y2 


The operation (8) causes 


pgs. 1 = Fs. (10) 
DiS spinors Y transform quite analogously under charge conjugation 
| yoy, vor (11) 


provided we use a representation where the matrices y, Ya Уз have real elements and 


ў 
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уа has imaginary elements. Gauge transformations, i. e. rotations in the plane (¢,, 65) 
through the angle х influence bosons and Dirac particles in the same way. 
ф-> eit Ф, р = ае y. (12) 
The fact that the transformation laws (10), (11), and (12) are of the same form for 
charged pions and for Dirac particles indicates that both are vectors in the t-space, 
i. e. both possess isospins 1. The only difference is that in the case of pion its isospin 
has three orientations —1, 0, 1 whereas in the case of Dirac particles there are only 
two orientations —1, 1 parallel to the third axis. Therefore pion possesses isospin 
of the type / but Dirac particle possesses an isospin of the type K. This is connected 
wth the fact that Dirac particles are subject to a Zitterbewegung but pions are not. 
Dirac spinors are spinors in x-space but vectors in £-space. On the other hand, 
Majorana particles) are spinors in x-space but scalars in t-space. 


Conclusions and summary 


The above considerations provide a tentative explanation, from a uniform 
viewpoint, of the problems listed in the introduction. An extreme usefulness of the 
conception of isospin implies that it cannot be merely a formal shorthand but is to 
be considered so seriously as the ordinary spin. This implies that an isospace possesses 
a physical reality similar to the reality of the usual space. It is possible to extend the 
law of conservation of the third component of angular isomomentum for weak inter- 
actions at the cost of introducing a composite, atom-like structure of particles, i. e. 


to characterize them by two different isospins J and K. The characterization of particles 
according to the quantum numbers b, K?, Kg, ГР, I, leads to a table of particles with 
striking regularities, This table suggests a close connection of the problem of rest 
masses with that of rotation in isospace. The particles resemble the neutrino of Lee 


and Yang in so far as their isospin K is necessarily oriented parallel (antiparallel) to 
a privileged direction for particles (antiparticles). This peculiar property implies 
isospace to be a euclidean subspace of a more-dimensional ps-euclidean space. It is 
assumed that both, ordinary space and isospace are two three-dimensional subspaces 
of a more-dimensional ps-euclidean space whereby the time axis plays the role of the 
privileged third axis in isospace. This explains why by a reversal of time the roles of 
particles and antiparticles can appear interchanged. Our assumption enables to under- 
stand why pion is ps-scalar. The reason of a restriction of orientations of the spin K 
can be traced back to a Zitterbewegung. 

The possibility of explaining why pion is ps-scalar seems to us very promising. 
We hope that the amalgamation of the *plane of the electric charge* with space-time 
into a more-dimensional world will enable to solve several further problems, in par- 
ticular to fulfil the old program of “geometrization“ of the electromagnetic field. 

The discovery of remarkable regularities among elementary particles seems 
promising. However, there are several unsolved questions. In particular, the table 2 
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suggests the existence of further unknown particles. Do they exist? There must be 
some selection rules prohibiting the existence of most of them: we have to avoid 
particles with half-integer charge as well as with charge larger than unity. There 
remain the following possibilities: a baryon with I = К = 0, a baryon with I= 0, 
К = 2 and a heavy meson with J = 0, К = 1. The baryon with I = К = 0 cannot 
exist since it would possess a small mass and all baryons would decay into it. On the 
other hand an existence of a baryon with J = 0, K = 2 is not apriori excluded. It is 
a relative of the hyperon A with a mass about 2860 electron masses. It would be nega- 
tively charged and could decay into the hyperon & with emission of a neutral pion. 
The heavy meson with Г = 0 and К = 1 should be negatively charged (but it should 
possess a positively charged antiparticle) and probably is essentially heavier than 
the usual heavy meson K. If its mass were sufficiently high, it could decay into a K- 
-meson with emission of a pion (or muon and neutrino). The existence of such a heavy 
meson was suggested by several experimental research workers. 


КРАТКОЕ СОДЕРЖАНИЕ 


E. Райски, О возможности связи между пространством и изопостранством. 


Сохранение третьей слагаемой изомомента (момент количества движения 
в пространстве изотопового спина) может быть обеспечено даже для слабых 
сопряжений за счёт предположения структуры частиц изопространства (два 
независимые изоспина I, а также К). Можно классифицировать частицы CO- 
гласно квантовым числам b, K2, Кз, 12, 13, где b означает число барионов. Эта 
классификация позволяет разместить частицы в таблицу, отличающуюся пра- 
вилами, которые нас интересуют в данном случае. Массы покоя можно объяс- 
нить с помощью этих квантовых чисел. Существует тонка и сверхтонка струк- 
тура массовых уровней определена 12, 13. 

Аналогия спина К (только две ориентировки) подсказывает мысль, что 
изопространство является эвклидовым подпространством некоторого псевдо- 
эвклидого пространства. Предполагается объединить изопространство с обыкно- 
венным пространством в псевдоэвклидовое пространство, которое имеет больше 
измерений. Тогда можно истолковать аномольный характер спина К с помощью 
,Zitterbewegung” (дрожащего движения). 

Предполагает 12-мерное пространство легко понять почему антипластицы 
можно считать, как частицы бегущие вспять по времени и почему отвес является 
псевдоскалярным. 
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ON THE INTERNAL DEGREES OF FREEDOM OF PARTICLES 
IN QUANTUM FIELD THEORY* 


By I. BIALYNICKI-BIRULA 


( Received Januar 16, 1958) 


In this paper we give the definition of the operators of the polarization and the 
magnetic moment in quantum field theory. For the free fields these operators are the 
constants of motion. We investigate exactly their eigenvalues for a spinor field. In the 
case of electrodynamics and mesodynamics we obtain the time dependence of these 
operators for the interacting fields. 


Introduction 


The purpose of the present paper is to investigate and to describe the internal 
degrees of freedom of particles in terms of the field operators of quantum field theory. 
Such operators have not been defined so far. They have to describe, among other 
things, the polarization and magnetic moment of the field. They must be the constants 
of motion in the case of free fields. In what follows we shall deal with particles having 
definite momenta. This means that the operators of creation and anihilation will appear 


as coefficients in the Fourier integral: 


v) = N f ar (u (да, (ре FF + Y (i oye 7) 


The momentum alone, however does not, in general, completely the state of the particle. 
The independent states having equal momenta are labeled by an index r in (1). 
In our paper we find an operator G, such that the one-particle states of the field: 
um О 000) |o > are eigenstates of С, belonging to different eigenvalues 
(gC? (p), respectively) for different values of the index r. Then the set of operators: 
the energy-momentum of field Р,, the number of particles № and the operator Go is 
the complete set of commuting observables for one-particle states of field. Having the 
description of the one-particle state of the field by means of eigenvalues of the operators 
Р, Мапа Go, we can describe all states of the field with the help sf the occupation 
ана This means that, Ше state vector of the field can be constructed as a linear 


superposition of the vectors: lo >, aj (р) |° >, b+ (р) [о >, af (pb; (р) |о > 


previous papers by the author (Bialynicki 1957 abc). 
(153) 


1 The present paper is based on 


у 
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1. The free fields 


1. Preliminary concepts and notation 


We shall investigate the fields, which satisfy the free field equations in the form (1). 
For the time being we restrict our consideration to the charged fields. 


(i h p" 3, + Mo) у (x) = 0 (1) 


(It is well known that every set of linear differential equations can be reduced to the 
form (1) by the introduction of the additional variables y, instead of derivatives). 
The properties of the field is greatly influenced by the form of the matrices В", as was 
pointed out by many authors (Umezawa 1956, Bhabha 1949, Harish-Chandra 1947). 
In the present paper the symmetrized and antisymmetrized products play an important 
role. We introduce a special notation to describe them: 


баир — (k гуа 21 p^: 2 per (2) 
perm 
Ant — (kN) DA p^ ... pur (3) 
perm 
k(k — 1) 


У) means the sum over all permutations of the indices. We have 5 = 1 when 
perm 


: k(k — 1 
is odd and S — 0 when ( ) is even. It is worthwhile noting that 48-6 = 0 


when k > 4. In what follows we shall assume that the field describes particles with 
definite masses M. From Eq. (1) the Klein-Gordon equation must follow for all compo- 
nents y,. This means that there exists a differential linear operator A(ihd) such that 


A (9) (ihB* 3 , + Me) = 129 де + M*c? (4) 

‚Тһе operator A can be constructed (Umezawa 1956) from matrices В“. It has the form 
A (ih8) = (ih)" (Mc udo. e DRE +... Ба", +aMc, (5) 

where «= 1, at = (--)"(К1) 1 > [Bet па Bok == йл ieee ase] (6) 


perm 


The equations (1) follow from the Lagrangian 


Е E 
L= = thpped wp + Мору, (7) 
pony T po, т) 9,9 Y y — yq. 
It is well known (Umezawa 1956, Bhabha 1949) that В“ and 7 must fulfil the following 


conditions on account of their hermicity and an invariance with respect to the Lorentz 
transformation: 


ПЕ орар: (8) 
A BRA = аў BY. Ау as 
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Л is the matrix transformin 

7 ) S g y to the new system of coordi ; this 

ЕО. Y y dinates; this means that 
= Ду x" = а ,х. Owing to condition (8) the expressions ри #8 y and 

— tt... : و‎ 

pA“ ТЕ y are tensors. We shall show that they are hermitian. 


(58) = (k!) > Вт!“ ... ВРО = y Str, 
(AMHR) + у) m (кл ( =й” ye (—)? Вт“ D" 8 = 


k (k—1) 


О Б у Е (9) 


. Hence it follows immediately that (p SY y)" = ySy and (y 40? y)* = AP y. 
With the help of Eqs (5), (6) and (8) it can be verified that 


At (ih?) y = 7 A (— то) (10) 
and that (то) (ihB* 9, + Mc) = (ihß"2, + Mc) A(ih?) (11) 
Let us now investigate the solutions of Eq. (1) in the form of plane waves р(х) = 


= 00 ( ра Ше! itive < 1 
р , which correspond to Ше positive and negative energy (px = 
= рх“, Po > 0). u‘*) (p) are solutions of the equations 


(Вер, + Mc) uP (р) = 0, (Вер, — Mc) up) = 0. (12) 


We label the independent solutions with the index r. Thus we have и (р) and uC?'(p). 
The index r relates to the internal degrees of freedom. We decompose the plane waves 


of the field y(x) with the help of these solutions: * 


ус) =N f aru (р) a, (pe №" + w (р) БА (р) eh), (03) 


АГ = 32. is an invariant element of the hypersurface p, p" = мее. N= Олт) ; 


Po . 
ap) and b} (p) are operators of the annihilation of a particle and the creation of an 


antiparticle, respectively. In the following sections we shall give an interpretation 
of the index r on the basis of quantum field theory. The normalization of the solutions 


иб" (p) is given (Appendix 3). 


2. Operators describing internal degrees of freedom 


To characterize the internal degrees of freedom of particles in quantum field 


theory we use the following operators: | 


4 T 
G(r) = He oe (14) 
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Tis a matrix acting on the index a of y,. We obtain the conservation law for the operator 


С(Г) with the help of the Klein-Gordon equation for free fields. 


д ih pao eee 
ie Oan E и 
Dia 
ih ih Е ae 
reri эл — — PLA | PLE = 
zx 1 er Ен (а 
E il verses A 
All operators having the form б, = 5 y 9, y, u.a да" are, of course, also con- 


а 


stants of motion. It can be shown in the momentum representation that the opera- 
tors С, are a nonessential generalization of the operators (14). The following 
equality must hold in order that the operator G be hermitian: 


[tn = qr. (16) 


The matrices Sk and Afı-#r have this property, as was shown previously. Thus 
we can construct a class of hermitian constants of motion G(S^r-^k), G(AM#R). 
We introduce a special notation for particularly important matrices At": 


ие = Ast — — (fn f pe po) 
1 vio 1 vA 
ди == 31 Ene Año => 31 er ? B, B3 De, (17) 


1 1 
Bs 7 4! 620 A io yes 4! ene, В» Pabe- 


If Г fulfils the condition (16) then the conservation laws for the operators (Г) follow 
from the theorem of Noether (Noether 1918, Hill 1951) on account of the invariance 
of the Lagrangian under the infinitesimal transformation: 


ду = ióA (Mc) A (#19) Гу. (18) 


Strictly speaking, the transformation (18) changes the Lagrangian (7) by adding a four- 
-divergence. (App. 1). We now pass to the momentum representation. We can obtain 
with the help of the formula (13) (App. 2) 


GI) = f ar fuer (р) Ги (p) а} (p) a, (p) + 
+ eu” (р) u™™ (p) bF (p) b, (p). (19) 
e = 1 for fermions and = = —1 for bosons. We introduce the following notation: 


АСЕ» (p! A) = Ars (p) EJU (p) I ue» (p) (20) 


nun —M— — 2 
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dra АСУ? (p) are hermitian on account of the hermicity of the operators С(Г) 
(in what follows we shall always assume that Г satisfies (16)). Now let us suppose that 


we pass from the set of solutions u+" (p) to the set of solutions ’u(#”(p) with the help 
of the transformation formulae 


‘Ww (p) = U (p) CO (p) 2 
Then the matrices A‘*)(p) transform themselves in accordance with the formulae: 
AD (p) = CD (p) AD (p) CH (р), _ 02) 


Owing to the hermicity of matrices A+? we can always diagonalize them by a suitable 
choice of unitary matrices ССР. For the solutions u‘*) (p) diagonalizing matrices AS (р) 
the formula (19) can be rewritten in the form 


с(Г) = [аГ{ =" (p) № (р) +28 (р) NO), (23) 
where g+” (p) are the eigenvalues of matrices AS (p), 
N+ (р) = а (р) а, (р), №? (р) = 5; (р) 5, (р). (24) 


In the formula (24) we do not sum over the index г. 

Let us investigate one-particle states of the field a+ (р) | 0», 6; (р) |0». We 
assume that the solutions и)" (р) determining the operators a+ (p) and b} (p) are chosen 
in such a manner that С(Г) has the form (23). We then obtain: 


G(T) a; (р) |0 > = ei? (paz (р) |0 > 

G (I) b? (р) |0 > = es (р) b; (р) |0 > (25) 
These states are thus eigenstates of the operator С(Г) belonging to the eigenvalues 
gr (p) and eg (p). If all eigenvalues are different (for definite p) then the operator 


С(Г) forms a complete set of observables for particles with definite momenta. Adding 
to the operator С(Г) the operator of energy-momentum 


p,— f 4 (p, DNO (р) + Pp 23 N; (P) (26) 
and the operator of the number of particles: $ 


м= [farce o-XNOQO . (27) 


we obtain a complete set of commuting observables for one-particle states of the field. 
That these operators commute follows from the equations: {№ (p), N» (р’)] = 0. 
The commutation rules for the operators С(Г) are obteined in App. 4. They have 


the form 


pryeri [aupa arado Фа (p) + 
— (Aq? дїн — AG? AMY BF (р) b, (P)}- (29) 
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If the matrices A, and A, commute for all momenta then the operators С(Г,) and 
G(T) commute. This fact is evident from the previous considerations. If the matrices 
A, and A, commute, then it is possible to choose a matrix C diagonalizing simulta- 
neously the matrices A, and Ay. This leads, on account of (28), to the commutation 


of G(T) and С(Г.). 
3. The free spinor field 
For the spinor field we have В" = y". The field equations take the form: 
(ћу, + Mo) y = 0. (30) 


A(ihd?) =(—ih y" 8 „+ Mc). The matrices 1, y", ив, II", y. form a complete set. All 
matrices 5А for k > 1 reduce to the product of the matrix y^ and 2%. We take 
as the point of departure the spinors и (p) in the form 


Up) ка u (р) = T | dup u” 1 
N = 2 Мор + Mo), КЕТ = (31) 
These solutions correspond to the following representation of уй: 

1 1 =; 1 
pl, һи=| qb У Е С |j. e» 
| "hi =] 1 

The matrices 4C? (Г) now have the form 

OEY AT: 

АЧ (y^) = (Moa, AO (y^) = (Mo pr 
A (и) = —(Мо-1рьд + [Mc(p, + МО (pa) p = — AC? (m), 

A (3) = (Mo (e x p) = —4O (3), (33) 
AG (по) = (Mc) po, AC? (x°) = (Mo) po, 
AP (л) = — о — [Mc (py + Mo]? (po) p =A (лу; 

АЧ (уу) = 0, A (yg) =0, 


1 
(и), = э Cm i", yk = کی‎ 


We can infer the following from the formulae (33) the operator G(Mcy") = p" is an 
operator of the energy-momentum of the field, G(— 1) is an operator of the number 
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of particles, and G(e) is an operator of the charge of field. The particle with the mo- 
mentum р = 0 has the matrices 4® in the form 


A) (и) = Фо, 4» (») —0, AD (до) = 0, AY E EM (34) 


Я D 7 ces A , ‹ 37€ 


arbitrary antisymetric tensor, n „is an arbitrary vector. The eigenvalues of the matrix 


Alu? пав) and A(z” n,) are calculated in App. 5. They have the form 


туз = + (Mo)! y^ D MD ID ns пе” М?с?, (35) 
S12 = + (Moy! Y (n^ py? medi». n" M? с?. (36) 


In the rest system we obtain from these formulae 


та = ЧЕ |m, |, Sa == з | | (37) 


The vector m, is a vector E, п" in the rest system of the particle and п, is the spatial 
part of the vector n^ in this system. If we assume that пр, = 0 апа лр, = 0 and 


that the vector 7, and the tensor n^" are normalized (i.e. $ n, n^' = 1 and n, n^ = — 1), 
then we obtain from (35) and (36) 
mis = +1, = + 1 (38) 


We shall now investigate more carefully the properties of the operator of the projection 
of polarization in the direction n,. Using the formulae (29) and (33), we can obtain 
the following relation (App. 6): 


[С (i y y" n), С G ys y" тру] | #2 = 0 (39) 
if and only if 
[(р„ пи)? = na n De M? с?] [(p, n)? = n» nO» М? с?] РЕ 
= [(n® ре) (и pr) — nP n M ср (40) 
| v > is one particle state of the field with a momentum p. Eqs. (39) give the necessary 
and sufficient condition that a particle have definite projections of the polarization in 


the directions no» and no» In App. 7 we show that the particle: described by the 


spinor и) (p) has a definite projection of the polarization in the direction s, = 


= UP (р) i ys y, Ш (p) and that s, Sl DSP 30: (41) 
Let us put пф съз and an — n, in formula (40): 


Кр, ne)? — n, n^ М? с?] = М? с? (n, s")?, (42) 
Using this equation in the rest system, we find that 


п = & ر‎ (Mc) Вр, (43) 
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is a necessary and sufficient condition that a particle have a definite projection of the 
polarization in the direction n,. The sufficiency of this condition is seen from the follo- 
wing: vector s, diagonalizes the matrix 4” (i ys y” S,), which means that 


А" (i ys y" s,) = u (p) i ys у" s, Ш (р) = — (0,)", (44) 
and therefore if we replace s, by n, defined by (43) we obtain 
A" (i ys y" п) = — a (0,)". (45) 


From (45) it follows immediately that the projection of the polarization in the direction 
n,isa = n,s". If a = + В then n, n" = 0 and n,s" = + (Mc) тр“. The projection 
in the direction of the light cone is equal to the projection of the momentum in this 
direction (up to the factor Mc). Our results can be summarized in the following manner: 
the states described by the orthogonal set of the spinors u‘*” (p) are the simultaneous 
eigenstates of the operator of the projection of polarization in the set of directions 
n, = as, + (Mc)! Вр, belonging to the eigenvalues 4,5 (п) = 4 a = +s,n* = 
= + (Mc) ү (пер)? — пл! M?c?. The upper sign relates to the state yc (р). We can 
examine the antiparticles in an analogous manner. 

Let us now study the operator of the projection of magneto-electric moment 
С (12 - y" п) from the same point of view In App. 4 it is shown that the tensor 
ти, = ИС (p) i ypu (p) diagonalizes the matrix А” (i[2* y"' m,,) and that 


ur’ 
1 m, m^" = land m*”p, = 0. The most general tensor which diagonalizes this matrix is 
a 


Пи» a AM yy 2n В (Mc) и t, x D, t). 


where t, is an arbitrary vector. In this case the conclusions are the same as before. 
The one-particle states of the field a^ (p) | 0 > relating to Ше set of orthogonal solutions 
ul ** (p) are the eigenstates of the operator С (i/2 y*”n,,) belonging to the eigenvalues 


412 (п „,) = а = + п, т" = + (Mo) үл“ Ray Py р" + in, n^" M?c? 


и 


if 
п, =: т, Ер В (Mc) (р, t, Py t) 


4. The free pseudoscalar (scalar) field 
We take the matrices f^ in the form 


мира pes EDS Eu age E S 


gio: E EE a M ИЕ 
ES В £ | 
Pom || , Ai SE ita = 
n да: il 23 
1 

N E B 


(46) 


DI A P 
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The solutions of the equations (8^ р, + Mc) и (р) = 0 are the following: 


Mc Mc 
ib ipo 1 — ipe 
Xeon экс үк E | 
pi |, и (р) = —= Tip (47) 
2M Р 
ү с — ips V 2 Mc +ip, 
— Pls +ips 
—1 
—1 
For the matrix у = 1 we obtain 
1 
1 
ао ом Tuc C» O (48) 


Using the formulae (46), we find 
n*= 0, В; = 0, 


m E | i г вв elis u” == клеш то ү = E 
| . . . . | . L . . . 
Using (47), we obtain WE) ав y = 0. All operators С(Г) describing the internal 


degrees of freedom vanish in this case. 


‚5. The electromagnetic field 
Let us consider the charged vector field in order to find the constants of motion 
С(Г) for the electromagnetic field. For this case we can apply the entire general theory 
previously dedueed. In the following we shall assume that the operators G(T) will 
have such a form as is obtained by passing from the charged field with nonvanishing 
mass M to the electromagnetic field with the help of the relations: ALS xe AD >A, 
иС в > Sap The Proca equations for the vector field 


Из а ne HE um 


ur 


G*”, DELI C uA“ (50) 
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can be written (according to Duffin and Kemmer) in the form 


(ihp*9, + Мо) y = 0, 


where 


pde 


(51) 


The matrices °°, л“ and В, for the vector field are given in App. 8. We construct from 


them the operators С(Г) according to the general theory: 


h 
бш) хэй 


where 


С (Bs) = 


би» == ER См. #0123 — 


ди 


дай (4008, Ar s с" д, 5, 


do? (4:8; С — 


do^ (eS С „+ С 


+1. 


CE) ^ if do^ (А; дл б — "д, A,), 


б» A), 


* 


Bv A 


do^ (A*« д, AF — A*8 д) Aa + Се д, ст — С*ё д, Gre), 


), 


(53) 
(54) 
(55) 
(56) 


(57) 
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Hence for the electromagnetic field we obtain 


G(f^) = f d * 4,8, p", (58b) 
1 > 
C, (ие?) = SE a 94 AP , (58c) 
1 Mr 
Се?) == | an for 5, fr, (584) 
G 1 ur Fe 
(85) = 19 ао Eo Ты, (58e) 
1 при 
C=} fam a Sj. (58f) 


We had to split the operator G(u*^) into the sum of two operators, because С; and С» 
have different dimensions. 


II. The interacting fields 
5. The Heisenberg picture 


Let us consider first of all two cases having the greatest physical importance: 

a. An interaction of the spinor field with the electromagnetic field. 

b. The charge independent interaction (three-component theory) of the spinor 
field with the pseudoscalar field. ч 


а. 


In this case we define the operators С(Г) by the following formulae: 


e= [Гуа (59) 
y, y are the operators of the spinor field in the Heisenberg picture. 


e 


бир = (дл — Аду — (91 + ic Any Y1 6 = y (60) 


In the above formula we replace 9,, which appears in Eq. (14), by D, according to the 

general principles (Corson 1951, Pauli 1941) as a result the operators (59) are invariant 

with respect to the gauge transformation. The equations satisfied by the field y have 
the form. 

(iy^9, +% + сув A) y = Dy = 0. (61) 

Applying the operator 4 

| D, = — iy^8, + × — yA, (62) 
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to Eq. (61), we obtain 
D,D,y = С + Aj) (-19в — eA") + x? + = y" 2 р = 0 (63) 
From Eq. (63) we obtain immediately 
y | — ib, +04 ión — eds фе + 22 = 0. (64) 


We rewrite Eqs (63) and (64) in the following form for the purpose of our further 
calculations: 


(9,9 — 214,9" — ie А", y) y = (ean EY — 5 yu sl Y, (65) 


VO ди + 2iü" An + ieAr,,) = Y | СА AP — x? — 5. 22 күб) 


We now find the variation of the operator С(Г) with time. 


OC) _ 
dal) — 


y (WT Dip)? = 5 yI(9,9^ — 2149, — Ар + 


E 5 v(9,9^ + 219,4? + іе, д). (67) 
Of course, we cannot put 4%, и = 0 in Eqs (65), (66) and (67), because A, are opera- 
tors. Substituting (65) and (66) into (67), we obtain 


И (68) 


Tkus the change in the operator С(Г) caused by а deformation of the hypersurface 
с(х) in the neighbourhood of the point x is proportional to the strength of the electro- 
magnetic field at the point x. The operator density у [Гуф is the coefficient of 
proportionality. Substituting into (68) the matrices 1, y^, iy^", iy; y", ys for Г, we obtain 
the following laws for the variation of the operators with time: 


Se = — efe фур, (69b) 
97 = ef pig", (603). 
Gly) _ (у (69e) 


да) ^" 


Mq m 9h 


In what follows we shall call the operators С(Г) the dynamical operators. In the case of 
a noninteracting field they are the constants of motion. For interacting fields they are, 
in general, functions of the hypersurfaces о(х). The operators р(х) y(x) (I-an arbitrary 
matrix) appearing in the right-hand side of Eqs (69) we shall call the kinematical 
operators. They depend upon the operators y(x) at the point x. They are the local 
operators (observables). Thus, to each matrix Г there corresponds one dynamical 
operator and one kinematical operator. Using this nomenclature we can read Eqs (69) 
in the following manner (assuming that the electromagnetic field f,, is the external 
field): 

a. The dynamical operator of the number of fermions is a constant of motion. 

b. The change of the momentum of fermion field is proportional to the strength 
of the electromagnetic field with the current as a coefficient of proportionality. These 
equations are related, of course, to the Lorentz equations: 


dp” _ 
ds _ 


359 efr,v’ 


c. The change of the dynamcial magnetic moment is proportional to the strength 
of the electromagnetic field with the kinematical magnetic moment as a coefficient 
of proportionality. 

d. It is easily seen from these equations that the change in the four-vector of the 
dynamical polarization is perpendicular to the four-vector of the kinematical polari- 


zation. 
e. The dynamical mesonic charge is a constant of motion in the electromagnetic 


field. 
b. 
In the case of the spinor field y in interaction with the pseudoscalar field we shall 
define the operators С(Г) by the formulae (14): 


(Г) = 5 1 PI Spd’. (70) 
а Р 
The fields y and y satisfy the equations in the Heisenberg picture in the form 


Cs 
(iyPd, +e + gyi Py = Dy=0, в= fe" (71) 


Applying, the operator 


to Eq. (71) we obtain 
D,D,y = 0,9" + 2 + e + igysy^,0;,) = 0. (73) 


Ud 
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From Eq. (73) we obtain immediately 
909,9* +12 + 80° + igysy" tD) = 0. (74) 
The dependence of С(Г) upon the hypersurface o(x) is given by the formula 


óG(I' i tm Е 
a "m 7 (Ге. (75) 
Using Eqs (73) and (74), we obtain 
óG(I' 
ми == 5 Pill уѕу итар (76) 


These equations have the same form as in the case of the electromagnetic field. The 
change of the dynamical operator is equal to the product of the kinematical operator 
and the strength of the field. We can now substitute for the matrix Г in Eq. (76) 
the product of the Dirac matrices and the matrices 1 or t;. Here we restrict ourselves 
to the matrices: 1, т» y*,1y*”, iyzy" ys. 


2 бы ie 

P = geyijD; iysy" tjv, (77b) 
E a = — go. фут, се 
A = — 8(0; Upiygy"y). (77d) 
a oe, SAS 
we = — gD; py" typ. (тї) 


Eq. (77a) gives the conservation law of fermions. Eq (77b) is the law of the change of 
the isotropic spin for fermions. We obtain the following dependence of the third compo- 
nent of the isotropic spin (charge) upon o(x): 


ÒG (та) Др CM 
UT = — E YPLY5 Y" YN Фи — gyn i ys y^ yp Фр. (78) 


Here yp, yy stand for the operators of the proton and neutron field and Ф = Ф, — 1Ф,. 
Eq. (77c) is the law of the change of momentum of the fermion field. This change is 
parallel to the strength of the meson field. Now the dynamical mesonic charge of the 
nucleon field is not a constant of motion. (Cf. (69e)). The change of this operator is 
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proportional to the current. Eqs (68) and (76) can be rewritten in integral form. We 
obtain 


E E 
EG sed tA Gesn LE [Г y] v dix, (79a) 


С (о) — C (a) = 5 | ФУ, ys y" til v dax. (79b) 


01 


6. The interaction picture 


In the interaction picture we define the operators С(Г) by the formulae (14): 


G(r) = г Гега (80) 


y and y satisfy the equations for free fields. Of course С(Г) are the constants of 
motion. In the interaction picture, however the operators alone are not important, 
their expectation values <y[o] | G(T) | v[c]7 = G(T) play an important role. The 
latters depend, in general, upon о(х). In the operators С(Г) there appear the time deriva- 
tives of the field operators and thus С(Г) are not related to the operators С.Г) in 
the Heisenberg picture (70) by an unitary transformation. Also the equality: 


<yle] | Go (Г) | plo] > = <Yy | Gy) | Vg i (81) 
does not hold. Using the definition of G(T ) and the Tomanaga-Schwinger equation we 
obtain: 

д GT) < ма | He, C) | viv] > (82 

да(х) h 


We can easily calculate the commutator [H, С] with the help of the following formulae: 
iver 
[G(R v (91 == y (9, — OLY), 


[6 (Г). F (| = — FOT Gr3, +9. (83) 
Assuming that the Hamiltonian H(x) has a form 
H(x) = v(9AGQ)vQ). (84) 


we obtain: : 


[G (7), На! = EGG), v 62] A0 v + (х) (а) [С (Г), ve. (85) 


y 
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and hence: 
[Н (x) SI] = 5 ФО Г, +) AQ) yl) + 


+ 5 9940 (79, - ЭГ. (86) 


We can rewrite Eq. (88) in the following form on account of the fact that y and y 
satisfy the free Dirac equations: 


UL (2), ED = 5 HAM T) 2, (ад + 


Deed с 
+ а vr APR). (87) 
The author is greatly indebted to Professor L. Infeld for his kind interest in this 


work and to Professor J. Werle for interesting discussions as well to Z. Bialynicka for 
assistance and discussions. 


Appendix 1. 


ih — = 
= y VB^9,y + Mc yy 


(A, 1) 
Let us perform the transformation: 
y > p= y + 104 (Mo) A (ihd) Py; (A, 2) 
then y >'p = y — 194 (Mo yr A (—ihd) (A, 3) 
on account of the сай (16) and Eqs (5) and 6) 
13 ^r. Si nr 
L= 5 ih уред, "р +Mc'p'y = 
lev ; lis - 
= 5 VhB^9,-- Мо) — > 'Y (ih B49, — Мо) "~ = 
e ees ;h gu 1 do 
==5 Ph pI, + Mo) y — g VGhf^8,— Mo y + 
1 . —/. . : 
+ 5 181 (Mo)? [р (ih f^ 9, + Mo) ARO) Пр + (А, 4) 


—pTA(—ihó) ih 8^8, + Mo) y + 
— wh вид, — Mo) A(ihd) Py + STA (— #9) (ih pu, — Mo) y]. 
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== means an equality up to (02)?. 


à.£ ='L — L = i6A(Mo)? [p (ih B49, + Mo) A(ihd) Ty + 


+ PTA (—ihd) (ih B» 8, — Mo) y] + 


14: * m 
— 3 102 (Me) [VA (— #9) ih Вер + Pih BLA Gh 8) y, 


From Eqs (4) and (5) we obtain 


8.6 = 194 (Mo)? [y (h?8,9» + M3c?) Гр — (29,91 + Мс?) ply] + 


a * > 
— 5 iA (Mc) [PTA (— гад) ihf^y + ріћв" A (ih) Гу, 
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(A, 5) 


(A, 6) 


| ES a A 
= 104 (Mc)! LL — > wl A (—iho) Вир — 2 wv BA Л 
7 


= 104 (Mo)! LoT P р о, | = 24. Е 


Thus the change of the Lagrangian is a four - divergence. It is well known (Hill 1951) 


that in this case we have a conservation law: 
IL 
M. mim 
G ES i d | 


IL dy = ЩЕ. Ви 164 (Mc)? A (ih) Гә + 


Mr 5 РГА (—ih9) i84 (Me)! B" y. 


From (A, 7) to (A, 9) we obtain 
О", = (p 9 Г) „= 0 
Appendix 2. 
Substituting in Eq. (14) y(x) in the form (13) we obtain 


cayo Pon fa far, fan 


ger (py ak (pJ eA Lr (p ke 49%) гӯ. 


utr (aga (p) e. хор) B (pe) = 


BR ‚far аг, · 
= 507% И пад 


(A, 7) 


(A, 8) 


(A, 9) 
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3 P 2+. -2= 

< [uC? (ру) Гис? (ра) (Pan + Pin) ar (ру) а; (ро) е аг 
E i (P3—P1) х 

+ a" (py) Га (рз) (Pan — Ры) br (p1) bs (py) ей Se 


4 3 (pep) x 
4- uC?r (ра) Dus (ра) (рзи Бри) аў (pa) Б: (ps) ей 


e = 1 (pop) х 
+U (py) Fus (py (pau — Pu) br (p) а: (ра) е # (А, 1) 


Taking the plane hypersurface с and performing the integration over d4x we obtain 


G(T) = 1 far, far, [uw ( p) Ги Оз ( pa) (Pao - Pro) 2 (py) as (p) д (Py --ра) 


+U (py) Ги: (рз) (—Peo — Pao) br (P1) ds (рә) 8 (P1 — Ра) + 


2i 
Є = ЖЕ — Ро 
+ ur (py) Pus (рә) (—Рэо + Pro) ar (ру) Ь (ро) д (Py + Po) ей a 


A x0 
+ ur (p) Ги» (ра) (Poo — Pro) br (p) as (р) 9 (а HBC ^ ]— 


= / аг [uC?* (p) Гиз (p) a (р) as (p) — Up) Fw (р), (p) ds (р)] 
(А, 2) 


We shall understand all expressions фГф as Ше normal products in order to obtain 
<G>p = 0. Hence Eq. (А, 2) has the form 


G(T) = far (p Pu? (р) аў (p) ap) + 


+ и)" (p) P'uC?* (p) b; (p) 6, (p)]. (A, 3) 


Appendix 3. 
From Eqs (12) we obtain 


Du (р) =a! DD (o) (Az 1) 


A proof showing that u? (p) В° и (p) and uC? (p) B°u (p) have the same sign 
for fermion field and differ in sign for boson fields can be found in the paper of Bhabha 


(Bhabha 1949) Eq. (8) allows for an arbitrary choice of the sign of the matrix y, and 
therefore we shall adjust it in such a manner as to satisfy u? B%u > 0 for bosons and 
fermions. From Eq. (12) we obtain 


U (p) B" р, и) (p) = —Mc u (p) u” (p), 


u (p) B" p, и (p) = Mc u™ (p) ш (p). | (Аз 2) 
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These equations give 
ur) (0) 8° и (0) — —'uC? (0) uC? (0), 


uC?(0)f* (0) = и“ (0) u (0) (Аз 3) 


in the rest system of the particle. 
Hence we have immediately: 


ut) (0) uC? (0) <0 
и) (0) и (0) |> 0 for fermions, (А, 4) 
< 0 for bosons. 


Eq. (A3 4) holds for an arbitrary p on account of the invariance of the expressions 
uS*) . u(9 , Thus we can assume that the solutions u( (p) are normalized in the follo- 
wing manner: 


иИО pju (p) = — б„ (Аз 5) 
u (р)шС (р) = вб„. 


Appendix 4. 
In order to calculate the commutation relations it is easy to use Eq. (19): 
(Г), 603] = f аг, / аг, А (АЯ) 
[u+ p) Ги (ра; (p3)a to) + 
+ вар) u py (Pb, ра), UO” (pT UA (ра)а (py) а„(рә) + 
+ eu ™™ (py Pou “(pabi (рә)®„(рә)] = 
= f dD dT, {up u (pa; (palp), и (pg) ри pas (рз)а,(ра)| + 
+ Hp) APA Y (ра), ИЗ YA YO (093 
We assume that these rules have the form | 
[а‚(}), a (Pl = PoP — р')д,,, (А, 2) 
[6,(Р), oF (PN. = pop — Рд, 


Hence 


[a (p) a (1), а  (ра)а (ра)! = (А, 3) 


= fax (pa, pa), (pr = p») di (paa (рз) д,„ӧ(Р1 — p2)| Po 
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and 
[b; (p))b,(p3). 0: (pb (ро)! = 
= [b} (o) ра), (р, — Pa) — bi (p2)b,(P1)5,p5(Pr — P2] Ро 
Substituting (А, 3) and (А, 4) into (A, 1), we obtain 
с), СОГ) = f ar - (atas m a (ар) + 
— (ASA — AD AYE (р) (р)}, 
where ИСЮ = 109 (pap), і = 1,2 
Appendix 5. 

Using Eq. (33), we get 

A (uns) = + (Mo) (ро + Me) (pn) (pa) = ра(оп) +a(p х п)) 
where (п) = пы» Nm = По, 
Ие matrix те obtain 

A? = (Ме)? [(Мс + p) °р? (p п)? + рат? + (рх n + 


— 2 (Mc — ро) po(P n) + 2pop (n x n) = 


= (Мо)? E nanag М?с? + n^? 77 = 12, 


(A, 4) 


(A, 5) 


(А, 1) 


(As 2) 


The square of the matrix A is thus the unit matrix multiplied by the number. T, 4* = 0, 


and therefore we get the following eigenvalues of the matrix A(+): 
Аз = + (Me)-! VV, n5 n*9 Mic? + n^ ni Pup 
For the matric 49 (n#p,) = 
= (Mà? [pang — Mean — (po + Mo)! (po) (рп) 
we obtain analoguosly: 
А? = (Mc)? [— M*?n,n" + (n^ p,)?], 
Hence: 


Aye = + (Mo)! V(n*p,)? — M*c*n, nt. 


Appendix 6. 


(As 3) 
(As 4) 
(A; 5) 


(А, 6) 


It is seen from (28) that |С(П.), С(Г.]а#(р)| 0 > = 0 if [A (р), AG? (p)] = 0. 
We shall show that this condition gives Eq. (40) when Г; = iysy"n®, Р, = iysy*n®. 
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Substituting the matrices AC? and AS” as in App. 5, we obtain 
[4GP, ASP] = [4,5, 4,5] = 2i0 (4, x А) = Ave. 
This commutator vanishes when с = 0, then also c? = 0. 
c? = (A, x Ag? = A? A? — (А, AS). 
А? апа А? are calculated in App. 5. Лов can readily be calculated: 
А, А, = — Mte nD ба + (p, ns) (pnb). 
Finally we obtain: 


а = Mey (руповв — nf nC? Мас] [(p,n 0»)? a Маса] — 
— (Mo)? [(n® p^) (nP р") — при M?e2]. 


Appendix 7. 


Let us consider two orthonormalized solutions of the Dirac equation и (р). 
We denote the vector И! (р) ѓууу^ uP! (p) by s". The traces of all matrices 
u C? (р) iygy^ u (p) vanish and hence we have 


s, = uO" (p)iysy ui?! (p) = — (p)iysy, u (p) (А; 1) 

It follows from the Dirac equation that s“p, = 0. It can be shown by simple calculation 
that s, 5^ = — (uC? (p) €?! (p)? = —1. Now let us construct the spinors 

uw (р) = iygy" s, Ш" (p) (A; 2) 


They satisfy the Dirac equation, because 
y^p wep) = —iygy pays Up) = —iysy's, Mou (р) = — Ме и” (р). (Ar 3) 
и 
; ГЫ : 
T (р) w (р) = iu? (p) ysy*s, uf? (р) = Ben araz 19 
WH (p) w (p) = WH (p) Гуру S175" За uH! (р) = 0,5,54 = да: (A75) 


From Eq. (A; 3) it follows that wt (p) = «и (р). From (A; 4) and (A, 5) we have 


= ай = btai 
al, = 1, a% = —1 and а 0, = д; and inte al, = ай = 0. Finally, we obtain 


wt (р) = uP (p), (р) = uC (р). (A, б) 


From Eqs (Ay 6) and (A, 2) we obtain: 
uC*(p)iygy"s,u Pp) =". (A, 7) 
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We shall consider in the magnetic moment an analogous way. Let us from the tensor 
M yy = И (pJiy „u Кр) = —u PA pyiy y up). From the Dirac equation it follows 
that m^" = p, = 0. It can be shown by calculation (de Broglie 1934) that 3 m,,, m*^ = 


= (uw? (р) vC?! (p)? = 1. 


Let +” (р) = 5 iy^" my, Ш)" (р) - 


1 
ураш) (р) = y iy! py y" тьш (p) = 
= — 5 Mey" mp, y! wr (р) = — Меш" (p). (A, 8) 
= т 5 АЕ 
и (р) wer) (p) = werr (p) wer? (p) ља ó,,. 


=== 


Finally, we have u” (p) $ iy" туи (p) = ( „)". 


Appendix 8 
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Appendix 9 


It is known that the operators of the free spinor field fulfil the rules of anticommu- 


tation in the form 


(v, (x), v (y)} = isa (x — у), (A, 1) 
where 56 (x — y) = (y^ 8,— x) Ах — y). (Ag 2) 
The function 4 is normalized in such a manner that 9, 4(x) | = — (x) 
From Eq. (A 1) we obtain 3079 


[C (7). v9] = 5 | f $ (=) P; y С) de^, y (x)] = 


, 


=. 14 
=> | s (v — 2) ГӘ y (и) da^ = — 5 (iy 9, — Ty (0. (463) 
с 


It is worthwhile noting that Eq. (A, 3) is consistent with Eq. (18). Taking the complex 
conjugate of Eq. (A, 3) we obtain 


[E 


(Г), BOM = — у PA Г@у*д„ +»). (Ao 4) 


КРАТКОЕ СОДЕРЖАНИЕ 


И. Биалыницкий-Бируля, О внутренных степениях свободы в квантовой 
теории поля. 


В этой работе мы даем формулировку операторов поляризации и магнитного 
момента в квантовой теории поля. Для свободных полей эти операторы не за- 
висят от времени. Исследовано собственные значения операторов поляризации 
и магнитного момента для спинорного поля: Полученная зависимость этих опе- 


раторов от времени для взаимодействующих полей в электродинамике и мезо- 
динамике. 


LITERATURE 


Bhabha, H., Rev. Mod. Phys., 21, 451 (1949). 

Bialynicki-Birula, 1., Bull. Acad. Polon. Sci. Cl. Ш, 5, 905 (1957a). 
Bialynicki-Birula, [., Bull. Acad. Polon. Sci. Cl. Ш, 
Bialynicki-Birula, I. Bull. Acad. Polon. Sci. Cl. Ш, 

de Broglie, L., L’electron magnetique, Paris (1934). . 

Corson, E., Tensors, Spinors nad Relativistic Wave Equations, London (1951). 
Harish-Chaudra, Phys. Rev., 84, 744 (1947a). 

Harish-Chaudra, Proc. Roy. Soc., 189, 372 (1947b). 

Hill, E. L., Rev. Mod. Phys., 23, 253 (1951). 

Noether, E., Gottingen Nachr., 235 (1918). 

Pauli, W., Rev. Mod. Phys., 13, 203 (1941). 

Umezawa, H., Quantum Field Theory, Amsterdam (1956). 


Vol. ХУП (1958) ACTA PHYSICA POLONICA Fasc. 2-3 


MULTIPLE PRODUCTION OF NONSTABLE PARTICLES IN PION 
NUCLEON COLLISIONS 


Ву V. $. BARASHENKOv Амр V. M. MALTSEV 


Joint Institute for Nuclear Research 


(Received Januar 17, 1958) 


The statistical theory of multiple production of pions, nucleons, nonstable particles 
and antiparticles in (л №) — collisions is considered by using the method described in 
(1), (2). The deductions of the theory may be put into agreement with the experimental 
data if it is assumed that "strange" and "usual" particles are produced in different space 
volumes. 


Introduction 


In papers (1), (2) it was shown that in the case of (NN)-collisions the relative 
multiplicity which is predicted by the statistical theory of multiple production and 
the charge distribution of secondary particles agree well enough with experiment when 
we are dealing with pions and nucleons. But it is in sharp conflict with experiment in 
the case of strange particles. In order to bring the aforementioned into agreement 
with experiment it was suggested that introduce one more parameterradius of the 
space-region in which the “strange” particles are “crystallized“ be introduced. It should 
be noted that the model under consideration is essentially different from that of Lepore 
and Neumann, (3) where the boundary of the space region is assumed to be diffuse 
although identical for different kinds of particles. | 

In this paper the model (1), (2) is applied to the multiple production of particles 


in (aN) — collisions. 


2. Results of Fermi-Belenky Model 


In Table 1 the theoretical and experimental values of the ratio o,/o, of the pro- 
bability of strange particle production to the probability of inelastic pion and nucleon 
production under the assumption that both strange (SA 0) and usual (5 = 0) particles 
are produced in the same space volume with a radius equal to the Compton pion wave 


(177) 


£ 
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length!. As well as in case of (NN) — collisions the theoretical value of this ratio many 
times exceeds its experimental value. An analogous result may be obtained also for 
the case in which we assume that all these particles are produced in a volume with 
radius of the order of the K-meson Compton wave length. Let us consider the mecha- 
nism of the multiple production of strange particles in more detail. 


3. The Results of the “Compound Particle“ Model 


Because of strong the interaction in a pion-nucleon collision a *compound particle" 
is originally produced. The fact that the nucleon "crystallization" starts simultaneuosly 
with pions from a volume whose radius is ~ (h/m,,c) is also accounted for in this strong 
interaction. Therefore only one parameter — the volume of the “erystallization‘ 
region? is included into the Fermi statistical theory. Quite another situation exists for 
sirange particles. We must consider that the pion interaction with K-mesons is conside- 
rably smaller than with the nucleons (otherwise о 0 is much greater than the experi- 
mental value, as was shown above). Owing to this the “flying away“ of free K-mesons 
will start from a smaller region with radius of the order of the K-meson Compton 
wave length. Then in the statistical weight formula the Fermi space factor V}, under- 
goes а change for V, or Vg?. 

The first case corresponds to then Schwinger and Gell-Mann hypothesis on the 
global pion interaction with baryons (7), the second case will occur if the pion interac- 
tion with A—; Z—; € — particles is considerably less than with nucleons. Аз can be seen 
from Table 2 the calculation with the weight factor V, (л-р) collisions is found to be 
closer to experiment (10)4. Since most of the strange particles are produced near the 
threshold one would not expect good agreement with experiment from the statistical 
theory. However, the minimum of the ratio о, [со in the region ~ 1 BeV, which is 
often of the statistical character, is confirmed by the experiment (10). 

In Table 3 the results are given of the calculation of the relative probability of 
possible reactions in (л-р) and (л+-р) — collisions with the energy Е = 5 BeV 
for the case of weight factors Г, and V4 (W,% and W,% respectively). In Table 4 
the corresponding results are given for (л+ р) and (a-p) — collisions with the energy 
Е = 7 BeV. In both tables the probabilities W, are expressed in percentages. The 
calculations are made under the same assumptions and using the same method as in (2). 


1 [n the cross-section of elastic production dp the cases of "elastic production" 


(л №) are not taken into account: 


of only one pair 


` (Na): в, = Y o (N- пл) + Xo (N-2N + пл) 


п-т n=0 


2 This fact is also reflected in the anomolous core diffusion in the nucleon. 

3 We use the notation as in (1), (2). : 

* An anologous calculation for (p—p) collisions at 3 BeV gives 16% for V, 5.7% for Va, 0.27% 
Гог Из. The second event is also found to be the closest to experiment ~ 3%. 
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In conclusion, we wish to thank Professor D. I. Blokhintsev for discussions. We 


are also grateful to Duan-I-Shi, V. L. Evteev and G. N. Tentyukova for assistance 
in the numerical calculations. 


Table 1 


Experiment 
0.95 24% 
1.3 19% TE 
5.0 134% 
7.0 201% 
Table 2 


Table 3 
(at p) (17 р) 
Reaction : 

n 2% W% W% W. 396 
№: пл 1 0.614 0.681 0,722 0,818 

2 8.68 9.67 10.2 11.6 

3 32.0 35.7 30.7 34.7 

4 317 35.3 30.0 34.0 

5 11.1 12.4 10.0 11.3 

6 1.42 1.58 1.26 1.42 
Ak - пл 0 0 0 0.246 0.0127 

1 1.64 0.0820 1.94 0,0984 

2 1.41 0.0696 1.66 . 0,0835 

3. 0.531 0.0260 0.536 0.0267 

ET ee O Е аа 

Zkena 0 ‚ 0.284 0,0143 0,335 0.0172 

1 2.88 0.144 3.40 0.173 

2 2.81 0.139 2.84 0.143 

3 1.02 0.0499 0.987 0.0490 
МЕЕ - пл 0 0.0753 0.0839 0.153 0.174 

1 0.415 0.463 0.418 0.473 

2 0.253 0,282 0.252 0.286 

3 0.0223 0.0249 0.0209 0.0237 
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Continuation of Table 3* 


(nt p) (пр) 
n W% W% 14 2% W% 
2N + Nnn 0 1:37 1.52 2.77 3.15 
1 1.51 1.68 1.44 1.29 
E-2k-n7 0 0.0176 0.03670 0.0357 0.02139 
1 0.0281 0.02107 0.0332 0.02129 
2 0.02175 0.04666 0.02185 0.04716 
AAN. na 0 0 0 0.0706 0.02126 
A * 2kk na 0 0.03178 0.05610 0.03361 0.04126 
ET Fe Ser Shi a sm it bin: CRT ERES 
E.2k-k-nz 0 0.03360 0.04123 0.03424 0.04148 
Table 4** 
PERRERA ARAS о ee 
(at р) (17 p) 
Reaction 
n W. 2% W. 3% W% W3% 
Nena 1 0.144 0.172 0.170 0.209 
2 3.06 3.64. 3.62 4.43 
t3 16.2 19.3 15.7 19.2 
4 24.8 294 23.5 28.7 
5 13.8 16.4 12.5 15.2 
6 3.07 3.65 2.72 3.32 
1 0.363 0.432 0.316 0.386 
Ackena 0 0 0 0.0595 0.02334 
1 0.622 0.0332 0.734 0.0403 
2 1.49 0.0785 1.75 0.0954 
3 0.959 0.0506 0.980 0.0526 
4 0.253 0.0131 0.245 0.0133 
ХЕ. пл 0 0.0691 0.02373 0.0815 0.02453 
1 1.21 0.0648 1.43 0.0788 
Зи 3.25 0.172 3.29 0.178 
3 1.86 0.0972 1.80 0.0967 
4 0.253 0.0131 0.234 0.0125 
МЕЕ + пл 0 0 0 0.0276 0.0342 
1 0.305 0,362 0.321 0.393 
2 0.545 0.648 0.550 0.669 
3 0.195 0.232 0.182 0.223 
4 0.0101 0.0120 0.02924 0.0113 


* The quantity in brackets signifies the number of zeros for example: 0.023 = 0.003 
** The quantity in brakets signifies the number of zeros. 
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Continuation of Table 4 


Reaction Sun) (=~ p) 
2 W% W% АД 7% 
2NN «nx 0 1.47 1:75 2.98 3.67 
1 10.3 12.3 8.55 10.4 
2 8.81 10.5 9.91 12.1 
3 0.486 0.431 0.204 0.249 
Sk na 0 0.02990 0.03405 0.0201 0.03850 
1 0.0495 0.02202 0.0584 0.02245 
2 0.0239 0.03969 0.0253 0.02106 
A 2kkna 0 0.03239 0.05873 0.03485 0.04183 
0.03104 0.05378 0.03122 0.05459 
У.о. Епл 0 0.02377 0.04137 0.03444 0.04167 
1 0.03171 0.05622 0.03181 0.05677 
E. Tha Y, — omo ме ТФ A de MA alta RE EF e 
AAN- nz 0 0 0 1.05 0.0204 
1 0.844 0.0149 0.433 0.02719 
УХМ-пл 0 1.43 0.0267 1.69 0.0324 
1 06% 0.02251 0.662, 0.02295 
а. ХМ. пл 0 0.960 0.0179 1.13 ' 0.0217 
1 0.572 0.02974 0.675 0.0118 
1-ZN-na 0 0.960 0.0179 1.13 0.0217 
à 1 0.572 0.02974 0.675 0.0118 
И RL LG ATE i ib Coli E ER M 
2.24 0 0 0 0.0741 0.04494 
725 0 0.0192 0.04123 0.0226 0.04149 
TAX 0 0.0783 0.04502 0.0925 0.04610 
ат 0 0.130 0.02241 0.263 0.02506 
A бе 8л NP D 
ANNE 0 0.0139 0.03732 0.0281 0.02154 
4 
ani 0 0.02693 0.03366 0.0141 0.03768 
ENNk 0 0.0157 0.0%828 0.0185 0.02101 
БИ? 0 0.02784 0.03414 0.02926 0.03503 
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КРАТКОЕ СОДЕРЖАНИЕ 


В. С. Барашенков и В. М. Мальцев, Множественное рождение нестабильных 
частиц при столкновении пионов с нуклонами. 


Статистическая теория множественного рождения пионов, нуклонов, неста- 
бильных частиц и античастиц в (л №) — столкновениях рассматривается методом, 
описанным в работах (1), (2). Выводы теории можно согласовать с эксперимен- 
тальными данными, если предположить, что „странные” частицы и „обыкновен- 
ные” частицы рождаются в разных пространственных объёмах. 
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THRESHOLD POTENTIAL FOR THE LOW FREQUENCY 
ELECTRIC DISCHARGE IN AIR USING A MAZE COUNTER 


Ву Н. D. SHARMA 
Department of Physics, University of Saugar, Saugor (M. P.), India 
(Received February 10, 1958) 


The threshold potential V,, for the low frequency electric discharge using Maze 
counters has been determined for air over a range of pressure p from 1 to 300 mm Hg. 
The V,,,— p curves are not linear, being concave towards the pressure axis. For a given 
wire diameter, 2r, V,, is seen to increase as the diameter of the outer glass cylinder 
increases. A comparison is also made using the rX-rp curves where X is the field strength 
at the wire surface. While these curves appear to coincide for low rp values, for higher rp 
values there is difference, the values of rX for the largest tube being the least. 


Introduction 


Coaxial cylinder coronas using metal electrodes under D.C. excitation have been 
studied extensively(!) because of their significance in Geiger counters. Studies using 
externally coated glass electrodes are now being made since 1946 after the discovery 
of Maze(?) that such counters also can work effectively. These studies have revealed that 
for these (Maze) counters (a) the Geiger threshold is slightly higher(?) and (b) the plateau 
extends to over an appreciably larger potential range(*) as compared to metallic Geiger 
counters. Under uniform fields also, the breakdown potential using glass electrodes 
is larger than that using metal electrodes(5). However, the Paschen’s law is obeyed for 
glass electrodes also. It was of interest therefore to see if the known relation 


И = гор (Ет): Дер) | 


for metal electrodes would also hold good in low frequency silent electric discharge, 
using glass electrodes. Here y, is the breakdown or threshold potential, r and R are 
the radi of the inner and outer electrode respectively, X is the field at the surface of 
the inner electrode, p is the ‘gas pressure and f denotes a function. 


Experimental 
The general experimental arrangement and procedure were essentially the same 
as reported earlier(6). The air used was carefully filtered through glass wool and was 
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stored over phosphorus pentoxide. The discharge tubes were essentially Maze counters. 
The dimensions of the tubes are recorded in Table 1. The tubes were excited by an Н.Т. 
transformer using 50 cycles A. C. The current, in the plate circuit of a detector, 
(1H5GT), which was connected across a serial carbon resistance in the L.T. line, was 
measured by a sensitive moving coil galvanometer. Before taking any observations 
the discharge tubes were aged for at least 30 minutes at every pressure under discharge 
well above the threshold potential, Vm. 


Table 1 


Dimensions of the tubes 


Metal and diameter of Outer diameter of the Inner diameter of the 
Tube No. 


the wire. cm. 2r glass tube. cm. 2R glass tube. cm 
1 Cu. 0.05 2.63 2.30 
2 Си. 0.05 2.03 1.82 
3 Cu. 0.05 1.41 Пе 
4 Си. 0.05 0.94 0.71 
5 Си. 0.05 0.785 0.55 
6 Cu. 0.05 0.69 0.45 
Results 


Values of the threshold potential, Vm, for different pressures of air in the range 
from 1 to 300 mm. Hg. are shown graphically in Fig. 1 (A), only for six tubes out of 
12 tubes used. The nature of the curves is seen to be similar to that for metal electro- 
des(t), the curves being concave towards the pressure axis. It is interesting to note that 
an increase in the diameter of the outer glass cylinder, keeping the central wire constant, 
leads to an increase in V,, as is to be expected. Fig. 1 (В) shows the variation of rX 
with rp, using the outer diameters of the outer glass cylinders for the calculation of X. 
Calculations of X were also made taking the inner diameters of the glass cylinders; 
however, they are not shown here due to their similarity with the present data. It is 
seen that while the гХ — rp curves appear to coincide for low values of rp, there is 
a slight divergence for higher values of rp. It is of interest to note here that for a given 
гр, TX goes on decreasing as the diameter of the outer cylinder is increased. 


Discussion 


The important role of the threshold potential, Vm, in reactions under electrical 
discharge and also in Ai phenomenon has been emphasized by Joshi(?). The threshold 
` potential, И „ is indicated by a sudden rise in the current flowing through the system, 
the initiation of a glow characteristic of the gas and the appearance of pulses on the 
current oscillogram. It is to be noted that an increase in the diameter of the outer 
cylinder diminishes the field at the central wire for the same applied potential, that is, 
to get the same field, a higher potential is necessary. An increase in the V m With an 
increase in the outer diameter of the tube is therefore to be anticipated. The physical 
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interpretation of this is simple. An increase in the cylinder diameter R tends to increase 


the effective gap length and hence tends to increase V. The results reported here 
are in accord with the view. 
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Fig. 1. (А) Variation of the threshold potential, V,,, with the pressure of the gas. The numbers on curves 
correspond to the numbers of tubes in Table 1. 


(В) гХ —rp curves for the same tubes. 


A comparison of data, using:a large number of tubes, as in the present case, may 
be made by using the relation $ 


pe Lie = Br LN E SOP a 
fog (RH Te C deu ^C 


Here JA, Va, V; .... Ri, Rs, Rg ..., and Xy, Xs, Хз ..., are the values of К, Y, and X 
for tubes 1, 2,3..., respectively. It is found that the rX-rp curves coincide within expe- 
rimental error except slight divergence at higher rp. This method of comparison has 
been used by a number of workers(8) for D. C. excitation using metal electrodes and 
the validity of this method has been discussed by Huxley (°). Essentially similar results 
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are reported by Jones and Williams(!9) recently for air under high frequency excita- 
tion. It is of interest to note that rX goes on decreasing as the diameter of the outer 
cylinder is increased, that for the largest tube used being the smallest. It may be pointed 
out here that contrary to the findings of other workers, who observed better agreement 
for higher pressures using metal electrodes, the spread of values in the present case 
is found to be appreciably larger for higher rp values i.e., for higher pressures. The 
present values of rX for the same rp as compared with those of Townsend(*) and 
Jones(19) are found to be higher and may be attributed to a fall of potential on glass walls. 

My grateful thanks are due to Dr. D.P. Jatar for suggesting the problem, his 
keen interest, valuable guidance and helpful criticism throughout the progress of this 
work. I am also greatly indebted to Prof. S. S. Joshi, D. Sc. (London) F.R.I.C. (Lond) 
F.N.L, F.A.Sc. for financial help which made this work possible. 


КРАТКОЕ СОДЕРЖАНИЕ 


Г.Д. Шарма, Пороговой потенциал злектрического разряда низкой частоты в воз- 
духе с использованивм счётчика Мазе. 

Определено с помощью счётчиков Мазе пороговой потенцил Vm электриче- 
ского разряда низкой частоты в воздухе в пределе давлений p от 1 до 300 мм Hg. 
Крывые Vm—p He являются линейными à вогнутыми относительно оси давления 
Для провода o данном диаметре 2r, Vm растет с ростом внешнего диаметра стек- 
лянного цилиндра. Сделанно также сравнение используя кривые тХ—тр, где X 
напряженность поля на поверхности провода. В то время когда для низких зна- 
чений тр кривые эти почти совпадают, для больших значений тр имеется неко- 
торая разност, TX для самого большого цилиндра имеет минимальное значение. 
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The basic assumption is that our world is a four-dimensional manifold of infinitesi- 
mal surface elements. The (two independent) basic vectors ig, ig tangential to the surface 
element are independent of the basic vectors ij, la, ig, iy, in the Mirkowski space but, 
together with them, span a six-dimensional vector space. The vectors 14, ig, is are real, i4, 
is, ig are imaginary. Rotations in the (is, i,)-plane mean gauge transformations. The 
vector space spanned by the imaginary vectors 14, 15, ig can be interpreted as the isospace. 
The privileged "third" axis in it is nothing else but the time axis. 

Such an extension of the physical space makes possible a geometrical interpretation 
of gauge transformations, charge conjugation, and isotopic spin transformations. It enables 
one to amalgamate the electromagnetic and mesic fields with geometry in a similar way 
as the gravitational field has been amalgamated with geometry. The extension of physical 
dimensions provides us with an explanation why pion must be ps-scalar. It accounts for 
the difference between the charge and the baryon number, and enables one to understand 
why only baryons (but not leptons) constitute the source of nuclear forces. 


1. A programme of geometrization 


Two distinct types of physical laws can be distinguished in traditional physics: 
one of them is closely connected with geometry, the other reaches beyond geometry. 
The laws of conservation of energy-momentum and of angular momentum express 
nothing else but homogeneity and isotropy of space-time and therefore belong to 
geometry. Some other conservation laws, e.g. the conservation of charge or of the 
baryon number have not acquired yet any geometrical interpretation. In other words, 
some of the groups of transformations (translations, rotations, and inversions in 
space-time) possess a geometrical meaning but some others (gauge transformations, 
charge conjugations, etc...) possess none. This situation is unsatisfactory. A conside- 
rable progress would be achieved if we succeeded in interpreting geometrically every 
transformation of importance in physics. 

In the contemporary field theory the gravitational field plays a distinguished role: 
it is the only field fitting naturally into the geometrical (Riemannian) framework. The 
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other fields constitute foreign elements in space-time. Numerous efforts in the past 
to geometrize the electromagnefic field have not been very successful. Nevertheless, 
it is generally believed that a considerable progress would be achieved if we succeeded 
in geometrizing further physical fields. 

Our programme is to find a geometrical interpretation of electric charge, baryon 
number, gauge transformations, charge conjugation, and to find a geometrical interpre- 
tation of the electromagnetic and nuclear forces fields including the notion of isobaric 
spin (isospin). In order to fulfil such a programme we shall extend the notion of physi- 
cal space to six dimensions by assuming that our proper world (Minkowski space) 
is embedded in a six-dimensional space with signature +++ — — — whereby the 
first four dimensions constitute the Minkowski space and the last two can be 
interpreted as an "electric plane" where rotations mean gauge transformations. 


2. Gauge transformations and charge conjugation 


Charged bosons are usually described by non-hermitian field quantities y. Their 
hermitian parts Фар sr given by the ET 


p= El Фа) + ipa), gt =, (Фа) — ipo) (1) 


can be interpreted as components of a vector in a plane to be called electric plane. 
Gauge transformations of the first kind 
ф = eg, фї = е (2) 
are nothing else but rotations in this plane. Charge conjugation 
port peo (3) 
can be interpreted as an inversion of the second axis in the electric plane. The same 
geometrical interpretation of gauge transformations and charge conjugation is appli- 


cable to fermions. The usual four-component spinors y can be split also into an hermi- 
tian and an antihermitian part 


y= A Vas + ipa), y! are (va) — iyo) (4) 


where фа) and y are hermitian. However, we cannot interpret them as components 
of a vector in the electric plane without some precaution since the decomposition (4) 
is ambiguous due to the arbitrariness of choosing a representation for the Dirac matri- 
ces. As is well known, the charge conjugation transforms spinors in a particularly 
simple way 


poy, pty | (5) 


provided a representation used where the elements of the y* matrices are real for 
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= 1, 2, 3 and imaginary for k = 4 (Majorana representation). As (5) is of exactly 
the same form as (3) and gauge transformations for fermions are of exactly the same 
form as for bosons, Dirac spinors can be also interpreted as vectors in the electric 
plane. Their components are given by the decomposition (4) in the Majorana repre- 
sentation. 
We conclude this section by stating that gauge transformations and charge conju- 
gation possess a clear geometrical meaning of rotations and inversions respectively 
in a plane to be called “electric plane”. 


3. A six-dimensional space. Isospin 


Question arises whether there exists any relation between the ordinary four- 
-dimensional space and the two-dimensional electric plane? The invariance under gauge 
transformations of the second kind (to be interpreted as a set of different rotations in 
the electric plane for different points in space-time) suggests that the two spaces are 
somehow related. It is suggestive to try a fusion of them into a common six-dimensio- 
nal space. 

A stumbling block for such a fusion seems to be the fact that the two additional 
dimensions obviously play a different role as compared with the ordinary four dimen- 
sions: we are not able to move freely in the electric plane so that only rotations but not 
translations in the electric plane seem to make any sense in physics. However, a similar 
objection could be raised against a fusion of time with the ordinary space into a four- 
dimensional space. Obviously time plays a different role as compared with space-like 
coordinates. We are not able to choose arbitrarily our position on the time axis but 
are compelled to be situated just at the actual “now” which is distinguished in compa- 
rison with any other position on the time axis by the mere fact that it is "now". We are 
able to choose arbitrarily our position and to make arbitrary translations only in the 
space-like dimensions whereas we can choose our position and make arbitrary transla- 
tions on the time axis only in our imagination but not in reality. In spite of obviously 
different role played by the time, the formal procedure of fusing it with space-like 
dimensions into a four-dimensional manifold was so rich in physical consequences. 
This enrightens us to hope that also a fusion of space-time with the electric plane 
into a six-dimensional manifold may turn out successful in spite of a distinguished 
role played by the fifth and sixth dimensions. i 


possesses a signature +++ ———. 

Let us stress once more the obvious fact (usually disregarded by physisists) that 
the present (our actual now”) is distinguished since it and only it means our actual 
reality. We know something about the past only via our memory and about the future 
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only via imagination, but past and future do not belong to our actual reality. Question 
arises whether "now" is a point on the time axis? This may be doubted. “Now” seems 
to be rather a very small but finite element dt of the order 10-23 or 107% sec. The 
inclusion of two further time-like dimensions suggests to introduce a generalized “now” 
as a three-dimensional element да? da? dx® of ultra-microscopic dimensions. We 
cannot avoid that only a single ultra-microscopic element in the time-like subspace 
constitutes our actual reality. This reality is changing constantly. The actual reality 
is disactualising itself whereas another actual reality emerges from the future. We are 
compelled to travel in an a priori prescribed way in some direction in the time-like 
subspace and our memory retains the picture of a one-dimensional sequence of actual 
realities forming a “tube”. Thus, two things are distinguished in our three-dimensional 
time-like subspace: the actual “now” and the direction of the sequence of different 
"nows" forming a tube. The points within the tube either belong or belonged or will 
belong once to our actual reality. The points beyond the tube never belong to our 
"now", they are passed by in our dramatic journey. In view of such distinctions it 
is natural to assume a system of coordinates with one axis x* along the "tube" and 
two further axes x5, x9 perpendicular to the tube. In such a frame of reference our 
reality (actual, past, or future) is confined to the immediate vicinity of 


Е E (6) 


We are allowed, in principle, to perform arbitrary coordinate transformations but 
the laws of physics do not need to be invariant under such transformations owing to 
the distinguished role played by the “tube”. Physically meaningful are only such trans- 
formations which (if regarded as transformations of space itself) are transforming the 
reality (actual, past, or-future) onto itself, i.e. the transformations mapping the imme- 
diate vicinity of (6) upon itself. 

In other words, we assume that our (past, actual and future) reality is not a four- 
-dimensional manifold of points but a four-dimensional manifold of two-dimensional 
(electric) elements. Every point in space-time is surrounded by an ultra-microscopic 
vicinity spreading out in the fifth and sixth dimensions. A vicinity of a point (6) in 
the electric plane is distinguished in a similarly inexplicable but similarly true sense 
as “now” is distinguished from "past" and "future". The only difference is that "now" 
in the fourth dimension changes constantly whereas “now” in the fifth and sixth 
dimensions remains always Ше same. The laws of physics have to be covariant only 
under transformations (linear or not) transforming the vicinity of (6) onto itself. 
Among such transformations there are Lorentz transformations and gauge transforma- 
tions. Nevertheless it is plausible that some of the properties of matter can be inde- 
pengeni of the anisotropy of our reality (brought about by the existence of a privileged 

"tube") and can possess a higher degree of symmetry with respect to the embedding 
six-dimensional space. For example, some properties of matter can be symmetric 
under arbitrary rotations in the time-like subspace (4, 5, 6) in spite of the distinguished 
role played by x‘. 
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There are two commuting groups of three-dimensional rotations in Ше six-di- 
mensional space: these in the space-like and those in the time-like subspace. The for- 
mer being connected with the ordinary spin, it is suggestive to connect the latter with 
the isospin. This interpretation of isospin accounts immediately why one of its compo- 
nents is intimately connected with the electric charge (rotations in the electric plane!) 
and why the law of conservation of isospin is sometimes violated: the reason is that 
the electric plane (or the time axis) is distinguished and that arbitrary rotations in the 
subspace (4, 5, 6) carry our proper world (our reality) into another world whose points 
never form our “now”. | 


4. Why ріоп is ps-scalar? 


Let us investigate a spinless boson field with isospin equal to unity. It is a vector 
field in isospace and, if our interpretation of isospace proposed in the foregoing section 
is all right, its components Фар Фор Фа) Should transform under rotations in the sub- 
space (4, 5, 6) like dx, dxf, dx* respectively. But, фр being a spinless field, neither of 
its components can be influenced by proper Lorentz transformations. However, x4 
being the time coordinate, Фа) cannot be a scalar but seems to be the fourth component 
of a fourvector in Minkowski space. This contradicts our assumption. The way out 
of this dilemma is as follows: we must assume фа) to transform like a product dx! dx? 
dx? dx* whereby both requirements are satisfied: pg, transforms like a vector compo- 
nent 4х4 under rotations in isospace and like a scalar under rotations in Minkowski 
space. But then the neutral pion is not a scalar but a ps-scalar in Minkowski space. 

Question arises as to the scalar or ps-scalar character of the charged pions. This 
problem is empty since we can reverse the signs of pajand Фо) by means of a suitable 
gauge transformation: a rotation in the plane (5, 6) through the angle л. Therefore 
we may call charged pions scalar or ps-scalar without any danger of contradiction. 
The only component that matters in this respect is the component фз describing 
neutral pions and this was shown to be unavoidably ps-scalar. Thus, we have the follo- 
wing theorem: a spinless field existing in three states of the charge must be ps-scalar. 

The possibility of explaining why pion must be ps-scalar constitutes, in our 
opinion, a crucial test of the fruitfulness of our idea of a six-dimensional space with 
signature +++ ——— embedding our reality: a Minkowski space together with 
its immediate vicinity in the fifth and sixth dimensions. 


5. Time reversal 


Let us introduce the operator of covariant integration over a hypersurface [doy 
This operator has, by definition, the formal property of a fourvector for the whole 
Lorentz group including improper transformations. Let us define the charge as a cova- 
riant integral of the charge and current density over a space-like hypersurface 


Q = J dorir à (7) 
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With this definition Q is a scalar. In spite of the fact that the charge density reverses 
its sign under time reversal, the charge does not. This definition of the charge is con- 
sistent with our geometrical conceptions according to which the charge appears to 
be connected with representations of the group of rotations in the electric plane, i.e. to 
be. proportional to the component of isospin in the direction-of the time axis. Isospin 
being an axial vector, its components cannot reverse sign under inversions of the 
respective axes. 

Let us denote the inversions of the axes in our six-dimensional space by /, (и = 
= 1,...,6) and their products by Г, = I,,. The electric charge is not influenced 
by I, but, according to the section 2, its sign is reversed by /¿. The geometrical inter- 
pretation of charge conjugation C in the case of tensors and four-component spinors 
is simply 


СА (8) 


Hereby the sixth axis is by no means privileged as compared with the fifth axis because 
I, differs from J, merely by a gauge transformation: a rotation in the plane (5, 6) 
through the angle z. 

In the framework of the six-dimensional world time reversal is ambiguous. We 
can reverse time either by means of an inversion J, or by means of a rotation in the 
subspace (4, 5, 6) around the fifth axis through the angle л. The latter is identical 
with the product ДД. Let us denote time inversion supplemented by a transposition 
by I;. It does not reverse the sign of the charge (though it reverses the sign of the charge 
density) whereas I,J, reverses the sign of the charge and is identical with time reversal 
in the sense of Wigner Т. As pointed out by Rzewuski (1957) the CPT theorem is 
a purely geometrical bussiness within the four-dimensional world. Indeed, in our 


notation CPT is identical with Поза: 


6. A geometrization of electromagnetic field 


Let us consider the fundamental metric form 
4 
ds? = > (dx*)2 — (dx>)? — (dx$)? | А (9) 
kim 


where an imaginary time coordinate is used (instead of g,, = — 1) in order to simplify 
computations. The above form is physically meaningful in the vicinity of х5 = x8 = 0. 
Lorentz transformations as well as gauge transformations of the first kind mean rota- 
tions preserving 15 = х6 = 0 and do not influence the shape of (9) But gauge transfor- 
mations of the second kind (also preserving x5 = х8 = 0) are non-linear transformations 


АК = Ё, д5 = کے‎ cosf --xSsinf, = قو‎ sinf + x6 cosf (10) 
with f = Дай) and convert (9) into 


ds? = g,, dx" ах (11) 
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where 
42 
En = Oj, — A f Se Es = See = —1 
= 6 Hs ps : 
Ens = fe Ere — X fe 85 = 0 (12) 


with the abbreviations 


12 = (x8)? +, fp =O, f (13) 
The contravariant components are 
HR 2 = 1+ (шугу 
ee р t fi (14) 


From (12) and (14) we compute the three-indices-symbols (Christoffel symbols) but 
write down only those which are of interest for our purpose 


Dis = fy I$ = Se Tis = Тв = Гі = 0. (15) 


The covariant derivatives of a six-vector possessing only nonvanishing electric com- 
ponents ф?, pé are 


9 = 9,9 HSP P = 9,9 — Љ Ф". (16) 
Going over to the usual complex notation for electrically charged fields 
1 1 
= — 5 +z 6 , T == = 5 — i 6 17 
Ф ya” QUT ya? $°) (17) 
уе get y 
Фа = (8, — И, Ф, = (9, + ifj ф! (18) 


which yields a geometrical interpretation of the well known expressions occurring 
in gauge invariant theories. 
The criterion of flatness of our six-dimensional manifold is the commutability 


of covariant derivatives. From (15) we get 
Фи Фи = OK, — 9, fj) = 0. (19) 


The gradient field f, = 9,f represents fictitious electromagnetic potentials analogous 
to the fictitious gravitational potentials introduced by means of a transition to curvili- 
near coordinates in Minkowski space. 

In close analogy to the theory of gravitation we can introduce a generalized rie- 
mannian manifold by dropping the condition f, = 9, f (11) and (12) and introducing 
in this way the true electromagnetic potentials. In this case the Christoffel symbols 


become 
1 1 
Tis =f— TAS fife, Dis — fa&— 5% x8 fi fki» 
Tis = Га = Гы = 0 0 
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where f,; = 9; f, — 9, f; denote the components of the field strength. Since we are 
confined to stay for ever in the immediate vicinity of the point denóted by (6) we can 
drop the additional terms appearing in (20) and get again expressions of the same form 
as (15) — (18). However, our space is no more flat since 


Ф;ы — Фдь = Ин. (21) 


Thus, the existence of an electromagnetic field f}; introduces into our six-dimensional 
space a deviation from the euclidean geometry. The well known expressions 2, — if, 
and 9, + if, occuring in the theory of electrically charged fields acquire a geometrical 
interpretation: they are covariant derivatives of vectors situated in the electric plane. 
The role of electromagnetic field in the six-dimensional space is quite similar to the 
role of the gravitational field in the four-dimensional space. Both types of fields cause 
deviations from ps-euclidean geometry and both determine the curvature. The electro- 
magnetic field becomes an intrinsically geometrical conception in a similar sense 
as is the gravitational field. 
The first set of Maxwell equations 


9 Ла = (22) 
is to be postulated (ог to be derived from other considerations) similarly as Ше equations 
of the gravitational field have to be postulated within the framework of general rela- 


tivity. Hereby the charge density plays a similar role to the energy-momentum density 
in the case of the gravitational field. 


7. The nucleon 


Our aim is to give a geometrical interpretation of the nuclear forces fields similar 
to that of the gravitational and electromagnetic fields. An introductory step towards 
a subsequent geometrization of gravitation and electromagnetism consisted in noticing 
an arbitrariness in the mode of description: a covariance under a group of transforma- 
tions involving some arbitrary functions of space-time coordinates. These are either 
general coordinate transformations in Minkowski space (for the case of gravitation) 
or gauge transformations of the second kind (for the case of electromagnetism). In 
order to geometrize the meson field it seems indispensable to look once more for an 
arbitrariness of the mode of decscription, i.e. for a new group of transformations more 
or less similar to the gauge group. Since. gauge transformations concern, in the first 
instance, the sources of the electromagnetic field we should concentrate now our 
attention on the sourges of nuclear forces, ie. on the nucleonic field. 

First of all let us remark that it is possible to describe interaction-free fermions 
with mass M by means of a modified Dirac equation. 


(у + убЕМ)у = 0. (23) 


The spinor y is of the usual type (transforms in the well known way under space and 
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Hereby the equation (33) goes over into 
Ге, + [1234(eM +70) Y = 0 (36) 


where the components of the vector С are (xem, Bem, yem). The theory is invariant 
under transformations of the type (35). 

Such a freedom of choosing a representation is analogous to the freedom of 
performing gauge transformations of the first kind. The next step will be analogous 
to gauge transformations of the second kind: We are allowed to assume different 
rotations of the type (35) for different points in space-time. Hereby the coefficients 
a, В, y, and consequently the components of the vector € become functions of x*. 
Still, its length is restricted for every point x* by 


(€(x*))? = m2. (37) 


Now, the components of this constant vector field expressed in terms of a variable 
basis can be reinterpreted as a variable vector field in terms of a constant basis. Finally, 
we can drop also the restriction (37) whereby the fictitious field e(x*) turns into a physi- 
cally meaningful meson field possessing three components Фар Фар Фа) in isospace 
and a ps-scalar character in the ordinary space. The interaction-free equation (24) 
or (36) is converted into a new equation including the interaction with pions. 

The above procedure of introducing the pion field in interaction with (nucleonic) 
matter is so much similar to the procedure of introducing the gravitational field in 
interaction with matter or the electromagnetic field in interaction with charged matter 
that we call it a geometrization. In contradistinction to the former cases it is not an 
encroachment upon the geometry of the four- or six-dimensional space but upon 
the geometry of the space (x!, ..., x9, s) where s is a spin variable assuming a discrete 
set of values. 

The above method of introducing the meson field was necessarily connected with 
the assumption of a complicated structure of the field constituting its source. This 
explains why neither bosons nor fermions with a simpler structure (leptons) can be 


sources of nuclear forces. 


КРАТКОЕ СОДЕРЖАНИЕ 


Е. Райский, Попытка геометризации мезоэлектродинамики. 


Основным предположением работы является положение, ' что физическое 
пространство является четырехмерным множеством бесконечно малых элементов. 
Базовые векторы iş M ig касательные к поверхностному элементу, вместе с базо- 
выми векторами пространства Минковского ij i» ig 14, напрягают шестимерное 
векторное пространство. Векторы i, ig із вещественные, a i4 i5 ig комплексные. 
Обороты в плоскости 15 ig означают ,gauge" трансформации. Напряженное про- 
странство на комплексных векторах ц 15 ig может быть интерпретировано как 
изопространство. | 

Это обобщение понятия пространства позволяет на геометрическую интер- 
претировку узоспина и электрического заряда. Электромагнитные и мезонные 
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поля становятся частью геометрии аналогичным образом, как гравитационные 
поля являются чисто геометрической концепцией (в общей теории относитель- 


ности). : 
В рамках этой теории получается естественное истолкование факта, почему . 


мезон л псевдоскалярный и почему только барионы (а не лептоны) являются 
исчотником ядерных сил. - E 
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time inversions) provided e is a ps-scalar assuming a value 1 in a certain system of 
reference but changing its sign under inversions. We choose this type of Dirac equation 
for the interaction-free nucleons. In order to describe proton and neutron by a single 
equation we have to double the number of dimensions in the space of the spinor index 
and replace (23) by 


(Гк, + Гиз М) = 0 (24) 


where Г* are Dirac matrices with eight rows and columns. We choose the following 
representation 


k у" 
r= | y Дита 4, (25) 
whence 
pue34 — pipepepa — г | (26) 
Фр 
As is well known, there exist some more matrices with eight rows and columns satisfying 


the usual anticommutation relations. Since we are interested in a six-dimensional 
space, we shall need two more of them, e.g. 


НЕ. 


It is easy to check that (25) and (27) satisfy 
ТЕ = 20, пу =E g | (28) 


The matrices (27) are connected with the fitfh and sixth dimensions in a similar 
sense as the four Г” are connected with the ordinary four dimensions. In particular, 


the matrix 
MIN 2 (Г5Г% — rers = (pen) (29) 


(where 1 means a unit matrix with four rows and columns) is connected with rotations 
in the plane (5, 6). Since 99 commutes with Г” it can be regarded as a constant from 
the point of view of Minkowski space. This shows that there are two independent 
conservation laws for the quantities 


N= [dere R= [ао РТР (30) 


where N is the nucleon number. 
The electric charge of the nucleon is given by 


Q= с: (N — В). (31) 
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The two summands of the charge imitate third components of two kinds of half integer 
isospins. 4R is nothing else but the third component of isospin in the traditional 
meaning of this word. Indeed, from (29) it is seen that Г [5,6] can be identified with 
the matrix 19) in isospace. But the quantity 4N is also connected with the group of 
rotations in the electric plane as was shown in the section 2. 

It is remarkable that the two summands of the charge (31) behave oppositely 
under inversions in the electric plane. The properly antisymmetrized N changes its 
sign but R does not so that J, transforms proton into antineutron. Thus, the electro- 
magnetic interactions involving nucleons violate the conservation of parity in the 
electric plane (and therefore also in isospace). This resembles closely the violation 
of parity conservation in the usual space by weak interactions. The term $(1 + y5) 
involved in weak interactions resembles closely the term $ (1 4-799) involved in 
the definition of the nucleon charge. In the case of weak interactions the theory is 
invariant under CP instead of P whereas in the case of electromagnetic interactions 
involving nucleons the theory is invariant under SP(e) instead of P(e) where P(e) 
means parity in the electric plane and S means charge symmetry, i.e. the operation 
transforming neutron into proton and vice versa. 


8. А geometrization of nuclear forces 


Assuming a representation (25) for the Г* matrices the equation (24) splits into 
two independent equations for the neutron and the proton respectively 


(99, + y Mp, = 0, (9, + Yy’eM)y, = 0. (32) 


However, in consequence of the electromagnetic interaction the masses of neutron 
and proton are not equal so that (32) should be written more exactly as 


(19, + уеМ фр = 0, (19, + *eM,)y, = О. (32^) 
These last equations may be fused again into 


[T*9, + Геза M + T:gm)] Y = 0 (33) 


where 
1 : 1 | 
М = s (М, + Mp), m= 3 (M, — Mp). (34) 


Now, the matrix "5: can be interpreted as the “third” matrix 1% of Pauli in isospace. 


We are allowed to change the representation in the spin index without changing the 
form of (25) 


— 


19 > ат y pi + yr, a? + B? + y? = 1. (35) 
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ON THE ELECTRON DENSITY DISTRIBUTION IN MOLECULES 
Ву WLODZIMIERZ Koros 
(Received December 28, 1957) 


Institute of Physics, Polish Academy of Sciences, Warsaw and Department of Electrochemistry, Institute 
of Physical Chemistry, Polish Academy of Sciences, Warsaw 


The results obtained in the previous papers (Kolos 1957a,b) for the excitation 
energies of ethylene and benzene molecules, and the astonishing agreement of the 
evaluated one-centre integrals (Kolos 1957c) with the semi-empirical values suggest 
that the study of other properties of the correlated orbitals 


y; (1,2) = N; p; (1, 2) по (1) 


may be interesting. In this note we shall present some conclusions resulting from the 
study of the electron density distribution corresponding to the wave function (1), 
in which q, (1,2) is the usual product or sum of products of one-electron orbitals 
and по denotes the mutual electron separation. 


It is evident that in the case of the helium atom the wave function (1) must give 
wrong results. If г, = 0, then for the second electron we have a 2s orbital (Slater- 
form) and as a result we must get too high an energy. For the hydrogen molecule the 
wave function (1) also gives a wrong electron density distribution. For example, in 
the V state, which is ionic in character, |y; (1, 2)|? achieves its maximum value (along 
the internuclear axis) when one of the electrons is close to one of the protons and the 
other is outside the molecule at a distance of about la.u. from the same proton. Thus 
the function (1), when applied to the hydrogen molecule, removes the electrons from 
the protons too strongly, and that is why it gives much poorer results than a function 
in which the correlation factor has the form (1 + А по). 

Consider now two electrons on a 2p, atomic orbital (with effective nuclear charge 
Z = 3.25) and a density distribution along the z axis. In this case |p; (1,2)|? has four 
maxima for 21 = + 22 = +0.615 a.u. If, however, the wave function (1) is used to 
calculate the density distribution, then the гу; factor cancels the maxima corresponding 
to гү = гу, and the other two maxima are obtained for „= — وچ‎ = +09 a.u. 
j. e. for a greater distance between the electrons and the nucleus. Besides, the density 


(201) 


y 
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function |y, (1.2)|3 has also four very small maxima, about 7 times smaller than the 
previous ones, when both electrons are located on the same side of the atom. The 
co-ordinates of one of these maxima аге z, = 0.4 au., =, = l5 ал., and the other 
maxima relate to the opposite side of the nucleus or to пита. of the electrons, 
respectively. 

If the wave function (1) is applied to the л electrons of the ethylene molecule 
similar results are obtained. For example, the V state remains ionic in character, but 
if one of the л electrons is on one side of the molecular plane the other must be on the 
opposite side. An increase of the average distance between the л electrons and the 
molecular plane is another consequence of the correlation factor гуз. 

Both these results seem to be reasonable and both are, to a certain degree, analo- 
gous to the results of other modifications introduced by various authors in the molecular 
orbital theory. For example, Araki (1956) in the calculations of the interaction energy 
between electrons eliminated the short-range part of the Coulomb potential l/r,,. 
From a formal point of view this may be regarded as analogous to the use of a wave 
function which prevents both electrons from being close together. On the other hand, 
Ohno and Itoh (1955) have shown that the results of the molecular orbital calculations 
for ethylene can be improved by using smaller values of the effective nuclear charge 
in the ionic wave function. Consequently the z-electrons are shifted to a greater 
distance from the molecular plane. The inclusion of the 3px atomic orbitals in the 
usual molecular orbital method (Jacobs 1955) also has a similar influence on the 


— 
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LETTERS TO THE EDITOR 


ON THE INTERPRETATION OF SLATER'S APPROXIMATION 
OF ATOMIC ORBITALS 


Bv ALOJZY GOLEBIEWSKI AND KACPER ZALEWSKI 


( Received December 21, 1957) 


In the theory of chemical bonds one assumes in general, (following Slater (1930)), 
the following form for the eigenfunction of an electron in the atom 


pie Nri" aire eerie ce (он), (1) 


where n* in this formula denotes the effective quantum number, Z* is the effective 
nuclear charge, and Y,, (9, p) is the angular part of the hydrogen-like function. 
Some authors of basic monographs on quantum chemistry (Gombash 1953, Eyring, 
Walter, and Kimball 1954) give an erroneous (in our opinion) interpretation of this 
approximation which is moreover not in agreement with the views of Slater (1930)!. 
Thus e. g. Gombash (1953) writes: 

»...Eigenfunctions (1) are a solution to the problem of the one-electron atom 
with a potential energy of 


LL tn (nte 1) 
= — I Па ее 2 
г + 2 e ag r2 ( ) 
...For the electron energy we obtain, by (1) and (2), 
1 Ze е 
eR E ба. 3 
2 Gg (3) 


It may be shown that function (1) is an eigenfunction of problem (2) for ещеп- 
value (3) if and only if I = 0, i. e. for s electrons. Thus, knowing the eigenfunction 
and the eigenvalue we can uniquely determine the potential corresponding to them 
from the Schroedinger equation. From this there results for (1) and | (3) the following 
potential: 


MS n* (n* — 1) —U1 +1 | 
as = T a Rude e (4) 
This is in agreement with 5 if and only if / 0. 


1 The same mistake is contained in Preuss's Tables (1957). 


. (199) 
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In the cited paper by Slater (1930) the justification for approximation (1) is not 
the same. On the basis of the results of Zener and the method of Hartree he arrives 
at the conclusion that the influence of the nodes in the radial port of trial function is not 
significant. The density distribution for hydrogen-like function and its term of the 
highest power in г, i. e. (1), is, in general, nearly the same. Slater therefore keeps the 
most important term of the hydrogen-like function. Eq. (3) is then an approximate 
expression for the EU As follows from (4), this formula contains an approxi- 
mation exept where n* = l + 1. 
In the monograph Quantum Theory of Matter Slater (1953) goes somewhat 
further by applying the rules which he gave for Z* and n* to hydrogen-like functions. 
Eq. (3) is then an exact eigenvalue in this approximation. 
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LABORATORY EQUIPMENT AND TECHNIQUES 


A NEW METHOD FOR THE RAPID MEASUREMENT OF THE 
SHRINKAGE FACTOR OF A NUCLEAR EMULSION 


By J. Bentsz AND С. Borowczak 
Higher Pedagogical School in Katowice 
(Received December 3, 1957) 


A rapid method for the measurement of the shrinkage factor of a nuclear emulsion 
by means of a suitably constructed pyknometer has been developed. The results obtained 
agree within the limits of error with the results from other methods known thus far. 


The rapid and accurate determination of the shrinkage factor plays an important 
role in nuclear emulsion technique. One method permitting a rapid and accurate 
determination of the volume of a given body is weighing. The shrinkage factor €s* 
has been determined by this method with the help of suitably designed pyknometer. 
We denote by М the mass of the unprocessed plate, by Mj the mass of the processed 
plate, by М, the mass of the glass backing, by Ма the mass of the pyknometer with 
the unprocessed plate and a given liquid, by M; the mass of the pyknometer with the 
processed plate and the liquid, by М; the mass of the pyknometer with the glass 
backing and the liquid. The volume of the unprocessed emulsion is equal to the 
difference between the volume of the liquid in the pyknometer containing the glass 
backing alone and the volume of the liquid in the pyknometer' containing the unde- 


veloped plate. We thus have 


M, + M, — М, — M. 
y, unprocessed emulsion — E 1 ma (1) 


where d, is the density of the liquid at a given temperature. 

Similarly, the difference between the volume of the liquid in the pyknometer 
containing the glass backing alone and the volume of the liquid in the pyknometer 
containing the processed plate determines the volume of the processed emulsion. 


d (203) 
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Hence 


MEM, АМ, ш 


pre D 
обесва emulsion d (2) 
с 


From the definition of “5“ we have 

V unprocessed emulsion  — М, 5 M, er M, == М, 
V processed emulsion Mi "s М, =— M, — М, 
If the emulsion is not mounted on glass, then (3) can be written as follows: 

Mo + M¿— M, 

M, + М — M; 


where Mg is the mass of the pyknometer filled with the liquid. 


(3) 


N ===: 


(3a) 


= 


Experimental Procedure 


The pyknometer should be so designed that the plate can be placed inside it. 
Fig. 1 shows three basic types of pyknometers that were used. 

The upper surface of the stopper that closes the pyknometer should be accurately 
ground so that the opening in the stopper can be covered by a ground glass plate. 
This prevents, to a large extent, the evaporation of the liquid. A guide line should 
be marked on the pyknometer; this decreases the error introduced by the failure to 
set into place accurately the upper part of the pyknometer. It should be noted that 
the accuracy of the measurements depends, to a large extent, on the concentricity of 


Fig. 1. 


the pyknometer. The following liquids were used in the measurements: distilled 
water and carbon tetrachloride (Oliver 1954); the latter, as we know, does not wet 
the emulsion. The results obtained for “5“ by means of pyknometer “a“ of Fig. 1 


are given in Table 1. The measurements were made for a 100 и type M-2 emulsion 
at a relative humidity of 70%. 
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Table 1 
Method of Measurement Shrinkage Factor 
(relative humidity 70% Sis 
during exposure and observation) 
Direct measurement of emulsion thickness before and after processing 1.93 + 0.04 
From dip of tracks (Jenny, Hurliman 1951) 1.94 + 0.08 
Pyknometer “а“ filled with water 1.932 + 0.107 
Weighing and taking the humidity into account (Rotblat, Tai 1951) 1.94 
Weighing and taking into account the humidity and the coefficient given by Rotblat 


(Rotblat, Tai 1951) 2,03 { 


In order to decrease the measurement error pyknometers “b“ and “с“ were made 
with larger ground joints. The results obtained with these pyknometers are listed in 


Table 2. 


Table 2 


Relative humidity of 85% during exposure and observation. 


Pyknometer “b“ filled with Pyknometer “с“ filled with 
H,O ССІ, H,O CCl, 
Shrinkage factor 2.051 2.100 Shrinkage factor 2.090 2.083 
+0.095 +0.116 +0.051 +0.068 


In the error calculation, errors made during weighing and in setting the upper 
part of the pyknometer into place were taken into account. 

As may be seen from Table 2, water is more suitable for the measurement than 
carbon tetrachloride (which evaporates quickly). On the basis of the results listed in 
Tables 1 and 2 it is proposed that pyknometer “c“ be used. This pyknometer has 
only one ground glass joint. 
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